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ERROR DIFFUSION USING NEXT SCANLINE ERROR 
IMPULSE RESPONSE 
Copyright Notice 

This patent specification contains material that is subject to copyright protection. 
5 The copyright owner has no objection to the reproduction of this patent specification or 
related materials from associated patent office files for the purposes of review, but 
otherwise reserves all copyright whatsoever. 

Technical Field of the Invention 
The present invention relates to the field of digital image processing and more 
10 particularly to apparatus and method for digital halftoning continuous tone images. 

Background Art 

Digital images 

An image is typically represented digitally as a rectangular array of pixels with 
each pixel having one of a restricted set of legitimate pixel values. Digital images may be 

15 black and white in which case the restricted set of legitimate pixel values is used to 
encode an optical density or luminance score. Digital images may also be colour where 
the restricted set of pixel values may encode an optical density or intensity score for each 
of a number of colour channels - for example Cyan, Magenta, Yellow and Black (CMYK) 
or Red, Green and Blue (RGB). 

20 Digital images are common where the image value per pixel per colour channel 

is an S bit unsigned number - providing intensity values in the range 0 through 255. Such 
images are often called "continuous tone" images because of the reasonably large number 
(256) of legitimate intensity values. 
Digital halftoning 
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By contrast, digital images are often printed on devices which provide a more 
limited variation in intensity or colour representation per pixel. For example, many 
Bubble Jet printers only provide the ability to print or not print a dot of each of a Cyan, 
Magenta, Yellow and Black ink at each pixel position. 

5 In order to print a digital image on a printer with lower colour resolution than the 

digital image, it is necessary to use the original image to generate an image with the 
required lower colour resolution such that the generated image has a similar appearance to 
the original image. This process of generating a digital image of similar appearance 
where each pixel colour value is within a smaller set of legitimate pixel colour values, is 

10 known as digital halftoning. 

For ease of explanation, digital halftoning is hereafter described for the case 
where the input image is a single colour channel, 8 bit per pixel image and the halftoned 
image is a 1 bit per pixel (bi-level) image. The input image values are known as 
greyscale values or grey levels. Extensions of digital halftoning from this monochrome 

15 bi-level case to cases where the halftoned image pixels have more than 2 legitimate 
output values (multi-level halftoning) and extensions to digital halftoning of colour 
images can be performed. See for example, "Digital Halftoning" Ulichney R., MIT Press, 
1987, pp 340-342. 

Consider the case of a single colour channel, S bit per pixel, image which is 
20 halftoned to generate a 1 bit per pixel (bi-level) image which is printable on a black and 
white bubble jet printer, for which each pixel of the halftoned image has one of 2 
legitimate pixel values, a "no dot" value and a "dot" value. 

An image region in the halftoned image will print as a minimum optical density 
("fully white") region when all pixels of the region have the "no dot" value corresponding 
25 to non-placement of an ink dot; it will print as a maximum optical density ("fully black") 
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region when all pixels of the region have the "dot" value corresponding to placement of 
an ink dot; and it will print as an intermediate optical density (halftone) region when some 
of the pixels of the region have the "no dot" value and some of the pixels have the "dot" 
value. 

A highlight (near white) region in the halftoned image will have only a few of 
the pixels with the "dot" value - in these regions the "dot" halftone output value is the 
minority or exceptional result. A shadow (near black) region will have only a few of the 
pixels with the "no dot" value - in these regions the "no dot" halftone output value is the 
minority or exceptional result. 

The role of the halftoning process is to generate the printable image so that an 
appropriate number of ink dots will print in appropriate pixel positions so that there is a 
good match between the optical density of image regions in the original image and the 
average optical density of the matching image regions in the printed image. 
Error diffusion 

Error diffusion is a digital halftoning method originally developed by Floyd and 
Steinberg and described in the publication "An Adaptive Algorithm for Spatial 
Grayscale", Proceedings of the SID 17/2, pp 75-77 (1976). Error diffusion as developed 
by Floyd and Steinberg is hereafter described as "standard error diffusion". 

An overview of standard error diffusion is now provided. 

The standard error diffusion algorithm processes pixels line by line from the top 
of the image to the bottom of the image. Each line (or "scanline") is processed one pixel 
at a time from left to right. 

The standard error diffusion algorithm employs a pixel decision rule in which a 
modified input image pixel value is compared against a threshold value. If the pixel's 
modified input value is less than the threshold, then the pixel's halftone output value is 
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assigned to be the lower halftone output value; and if it is greater, then the pixel's 
halftone output value is assigned to be the higher halftone output value, 

Following determination of the pixel halftone output value, an error is 
determined for the pixel as the difference between the pixel's modified input value and 
the pixel's halftone output value. The error is distributed to neighbouring, as yet 
unprocessed, pixels according to a set of weighting coefficients. 

A pixel's modified input value is the sum of the pixel's input image value and a 
neighbourhood error value for the pixel. The neighbourhood error for a pixel is the sum 
of the errors distributed to that pixel from previously processed neighbouring pixels. 

The set of weighting coefficients is known as an error diffusion mask. Each 
weighting coefficient of the error diffusion mask is associated with a pixel offset. When a 
pixel is processed, the error distributed from the pixel to an as yet unprocessed pixel is the 
pixel error multiplied by the weight corresponding to the offset from the pixel to the 
unprocessed pixel. 

Note that the error at a pixel which is distributed to its neighbours can be 
considered as a sum of the neighbourhood error at the pixel and the pixel-only error being 
the pixel's input value less the pixel's halftone output value. 

The error diffusion mask described by Floyd and Steinberg is shown in Fig. 1 . 
The error at a pixel is distributed in the proportions 7/16, 1/16, 5/16 and 3/16 to 4 
neighbouring pixels which are as yet unprocessed as indicated in the diagram. Fig. 1 
shows an error diffusion mask 100 in which an error for a current pixel 108 is distributed 
in the proportions indicated to four neighbouring pixels eg. 106 in the fractional 
proportions shown. Previously processed pixels eg. 110 are shown using a shaded 
representation, and it is noted that the current pixel 108 lies on a current scanline 102. 
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Previously processed pixels lie above a bold line 104 in Fig. 1, and a scanline processing 
direction is depicted by a horizontal arrow 112. 

It is noted that the sum of the weighting coefficients is 1. As a result 100% of a 
pixel's error is transferred to its neighbouring pixels. If the sum of weighting coefficients 
were greater than 1, then error would be amplified and could build up without bound. If 
the sum of weighting coefficients were less than 1, then error would be reduced. By 
having the sum of weighting coefficients equal to unity, the average intensity of a region 
of the halftoned image tends to match the average intensity of that region in the input 
image which is a very desirable characteristic of a halftoning process. 

In implementations of standard error diffusion where multiplication of error by a 
weighting coefficient results in significant rounding errors, it is necessary to co-ordinate 
the calculation of error distributions to neighbouring pixels so that effectively 100% of 
current pixel error is transferred on. This can be achieved by determining error 
distributions to all but one of the neighbouring pixels by multiplication by a weight and 
determining the error distribution to the remaining neighbouring pixel by subtracting the 
other error distributions from the total error to be distributed. 

For improved execution speed, error distributions corresponding to a particular 
pixel error value are often determined in advance and retrieved from a look up table. 

In the Floyd and Steinberg error diffusion mask, the pixels which receive error 
distributions from a current pixel are on the current scanline and succeeding scanline 
only. An implementation of error diffusion with this mask requires the use of a single 
"line store" memory to store neighbourhood error values. The memory is referred to as a 
"line store" because it is required to store a neighbourhood error value for each pixel of a 
scanline. Many of the modifications to error diffusion which are referred to below are 
achieved at the cost of extra memory for an additional one or more line stores. 
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Error diffusion algorithms are used in the printing and display of digital images. 
Many modifications to the standard error diffusion algorithm have been developed. 
Worm artifacts 

The error diffusion algorithm suffers from the disadvantage that in image regions 
of very low or very high intensity, it generates a pattern of image values in the halftoned 
image which are poorly spread - the exceptional values are concentrated in wavy lines. 
These patterns can be very noticeable and distracting to a viewer of the image - they are 
often known as "worm" artifacts. 

Fig. 5 shows a section of a halftone output image generated by Floyd Steinberg 
error diffusion which shows worm artifacts. This halftone output was generated by bi- 
ievel halftoning of an 8 bit per pixel monochrome source image with a constant grey 
value of 253. Fig. 5 shows a halftone output image 500 which, as noted, shows worm 
artifacts, as illustrated by, for example, pairs of pixels 502, and a quadruplet of pixels 
504. 

Modifications to error diffusion have been developed to reduce worm artifacts. 
Prior art methods of reducing worm artifacts 

One method of reducing worm artifacts is by addition of some randomisation. 
The randomisation may be achieved by adding noise to the input image, by adding noise 
to the thresholds or by randomising the error distribution to neighbouring pixels. A large 
amount of noise or randomisation can be added to fully avoid worm artifacts; however, 
this also seriously degrades the halftoned image. 

Another method of reducing worm artifacts is to vary the direction in which 
scanlines are processed. By way of example, U.S. Patent No, 4,955,065 titled "System 
for Producing Dithered Images from Continuous-tone Image Data" to Ulichney discloses 



J 91 201 



581700US.doc 



error diffusion with perturbed weighting coefficients and bi-directional scanline 
processing. 

Another method of reducing worm artifacts is by use of larger error diffusion 
masks. For example, in "Error Diffusion Algorithm with Reduced Artifacts", Eschbach 
R„ Proceedings of the IS&T's 45th Annual Conference, May 10-15, 1992, either a large 
or a small error diffusion mask is used depending on the input image grey level. The 
large error diffusion mask is suited to use in image regions of very high or very low 
grayscale, reducing the worm artifacts in those regions. A disadvantage of this method is 
that large error diffusion masks which distribute error to pixels of more than 1 succeeding 
scanline require additional error line stores and associated processing. 

A further method of reducing worm artifacts by use of an "extended distribution 
set" style of error diffusion mask is described in U.S. Patent No. 5,353,127, titled 
"Method for Quantization Gray Level Pixel Data with Extended Distribution Set" to 
Shiau and Fan. This patent describes error diffusion masks including pixel positions of 
only the current and succeeding scanline with additional pixel positions on the next 
scanline to the left of (that is, behind) the current pixel. 

Fig. 2 shows an error diffusion mask of U.S. Patent No, 5,353,127. An error for 
a current pixel 208 is, in this case, distributed to five neighbouring pixels eg. 206, where 
each of these five pixels receives a fractional proportion of the error as indicated. 
Scanline processing is from left to right as depicted by an arrow 212, and previously 
processed pixels, eg. 210, are full shaded and lie above a bold line 204. The current pixel 
208 lines on a current scanline 202. While this method is successful at reducing worm 
artifacts and only requires a single error line store, worm artifacts are still evident for bi- 
level halftoning of 8 bit grey scale image data, as can be seen in Fig. 6. 
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Fig. 6 shows example halftone output generated by bi-level error diffusion, using 
the mask of fig. 2, for an 8 bit per pixel source image with constant grey value of 253, 
Fig. 6 shows an image 600 which exhibits worm artifacts as exemplified by reference 
numerals 602 and 604. 
Prior art methods of preventing artifacts 

While the above methods are successful in reducing worm artifacts, complete 
prevention of worm artifacts can also be achieved. 

One method of preventing worm artifacts is by modulating threshold values. 
U.S. Patent No. 5,535,019 titled "Error diffusion halftoning with homogeneous response 
in high/low intensity image regions" to Eschbach discloses a modification to error 
diffusion which adjusts the error diffusion threshold according to the halftone output and 
according to the input intensity using a threshold impulse function, for the purpose of 
preventing worm artifacts. 

Another modification to error diffusion which prevents worm artifacts by 
threshold modulation is described in "A serpentine error diffusion kernel with threshold 
modulation for homogeneous dot distribution", Hong D., Kim C, Japan Hardcopy '98 pp 
363-366, 1998 which is also published in IS&T's Recent Progress In Digital Halftoning II 
(1999) pp 306-309. 

Another method of preventing worm artifacts is the addition of grey level 
dependent periodic noise to the input image, described in "Modified error diffusion with 
smoothly dispersed dots in highlight and shadow", Japan Hardcopy '98 pp 379-382, 1998 
which is also published in IS&T's Recent Progress In Digital Halftoning II (1999) pp 
310-313. 

Another method of preventing worm artifacts is by imposing output position 
constraints. An example of this method is provided in "An error diffusion algorithm with 
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output constraints for homogeneous highlight and shadow dot distribution", Marcu G., 

Proceedings of SPIE, Vol 3300, pp 341-352 (1998). 
Disadvantages of nr ior art methods of nrevenrine worm artifacts 

In many cases the desirability of a halftoning algorithm is determined by how 
fast it executes and how easy it is to implement. For example, in software 
implementations in a printer driver on a general purpose computer, the algorithm 
execution speed is very important; whereas in special purpose hardware, the algorithm 
complexity and memory usage are very important because they relate strongly to the 
expense of the circuitry. 

Use of additional line store memory by a halftoning algorithm generally 
indicates that it will execute slower in software and is more expensive to implement in 
hardware. 

The modification to error diffusion disclosed in US Patent No. 5,535,019 to 
Eschbach requires use of memory for an additional line store to store threshold 
adjustment values generated by preceding scanlines for use by a current scanline. The 
modification also requires additional processing including addition of threshold impulse 
values and dampening of threshold values transferred to subsequent scanlines. 

The modification to error diffusion described in the previously mentioned paper 
titled "A serpentine error diffusion kernel with threshold modulation for homogeneous 
dot distribution" also requires use of an additional line store memory to store threshold 
adjustment values. Additional processing is also required to diffuse threshold adjustment 
values. 

The modification to error diffusion described in the previously mentioned paper 
titled "Modified error diffusion with smoothly dispersed dots in highlight and shadow" 
requires memory for additional line stores to store several neighbouring scanlines of 
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processed input image data from which filtered input image values for a current scanline 
are determined; the filtered input image values arc used in turn to determine the noise to 
be added to the input image data. 

The modification to error diffusion described in the previously mentioned paper 
titled "An error diffusion algorithm with output constraints for homogeneous highlight 
and shadow dot distribution" requires memory for additional lines stores to store 

halftoned image data for several previously processed scanlines. This modification also 

includes processing to exclude minority halftone output results when that result is present 

in a certain portion of the previously processed scanlines. 

In summary, all the modifications listed above which prevent worm artifacts in 

error diffusion require use of additional line store memory together with the processing 

associated with use of that additional memory. 

References made to prior art documents in the present description in no way 

constitutes an acknowledgment that the prior art documents are part of the common 

general knowledge. 

Summary of the Invention 

It is an object of the present invention to substantially overcome, or at least 
ameliorate, one or more disadvantages of existing arrangements. 

According to a first aspect of the invention, there is provided a method of halftoning 
an image, said method comprising steps of: 

determining an output value of a current pixel on a current scanline using a sum 
of an input value for the current pixel and a neighbourhood error value at the current 
pixel; 
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determining an error at the current pixel as the difference between (i) the sum of 
the input value for the current pixel and the neighbourhood error value at the current 
pixel, and (ii) the output value of the current pixel; and 

adding a proportion of the error at the current pixel to neighbourhood error 
5 values at as yet unprocessed pixels of a subsequent scanline in accordance with a next 
scanline error impulse response; wherein said next scanline error impulse response: 

approximates a function which spreads with self-convolution in proportion to a 
degree of self-convolution. 

According to another aspect of the invention, there is provided a method of 
10 generating an error diffusion mask suitable for use with any of the aforementioned 
methods. 

According to another aspect of the invention, there is provided an error diffusion 
mask suitable for use with any of the aforementioned methods. 

According to another aspect of the invention, there is provided a method of 
15 halftoning an image, said image comprising a plurality of pixels each having an input 
value and an assignable output value that can take on one of at least two output values, 
where pixels are processed scanline by scanline, and scaiilines are processed one at a time 
from the top of the image to the bottom of the image, and where a scanline is processed 
pixel by pixel either from left to right or from right to left, and where the processing for 
20 each pixel comprises the steps of: 

(a) determining the output value of a current pixel using a sum of the input 
value of the current pixel and a neighbourhood error value for the pixel; 

(b) determining an error at the current pixel as the difference between, 
firstly, the sum of the input value of the current pixel and the neighbourhood error value 

25 for the pixel, and secondly the output value of the pixel; 
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(c) adding proportions of the error at the current pixel to the neighbourhood 
error values of yet to be processed pixels of the current and next scanline; 

and where the said proportions of the error at a current pixel are designed so that 
the next scanline error impulse response, being that function which maps 

(A) from a horizontal pixel offset; 

(B) to the total proportion of the error at the current pixel added to the 
neighbourhood error of that pixel of the next scanline which is displaced by the horizontal 
pixel offset from the current pixel, following complete processing of the current scanline; 

approximates a function which spreads with self-convolution in proportion to the 
degree of self-convolution. 

According to another aspect of the invention, there is provided a method of 
halftoning an image, said image comprising a plurality of pixels each having an input 
value and an assignable output value that can take on one of at least two output values, 
where pixels are processed scanline by scanline and scanlines are processed one at a time 
from the top of the image to the bottom of the image, and where a scanline is processed 
pixel by pixel either from left to right or from right to left, and where the processing for 
each pixel comprises the steps of: 

(a) detennining the output value of a current pixel using a sum of the uiput 
value of the current pixel and a neighbourhood error value for the pixel; 

(b) determining an error at the current pixel as the difference between, 
firstly, the sum of the input value of the current pixel and the neighbourhood error value 
for the pixel, and secondly the output value of the pixel; 

(c) selecting, using the current pixel input value, a set of proportions and a 
set of corresponding pixel position offsets, from a family of sets of proportions and 
corresponding pixel position offsets; 
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(d) adding the selected proportions of the error at the current pixel to the 
neighbourhood error values of yet to be processed pixels at pixel positions offset from the 
current pixel by the selected corresponding pixel position offsets; 

and where each set of the said family of sets of proportions and corresponding 
5 pixel offsets, is designed so that the next scanline error impulse response corresponding to 
that set, being that function which maps 

(A) from a horizontal pixel offset; 

(B) to the proportion of the error at the current pixel added to the 
neighbourhood error of that pixel of the next scanline displaced by the horizontal pixel 

10 offset from the current pixel, following complete halftone processing of the current 
scanline using only the said set of proportions and corresponding pixel offsets; 

approximates a function which spreads with self-convolution in proportion to the 
degree of self-convolution. 

According to another aspect of the invention, there is provided a method of 

15 halftoning an image, said image comprising a plurality of pixels each having an input 
value and an assignable output value that can take on one of at least two output values, 
where pixels are processed scanline by scanline and scanlines are processed one at a time 
from the top of the image to the bottom of the image, and where a scanline is processed 
pixel by pixel either from left to right or from right to left, and where the processing for 

20 each pixel comprises the steps of: 

(a) determining the output value of a current pixel using a sum of the input 
value of the current pixel and a neighbourhood error value for the pixel; 

(b) determining an error at the current pixel as the difference between, 
firstly, the sum of the input value of the current pixel and the neighbourhood error value 

25 for the pixel, and secondly the output value of the pixel; 
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(c) selecting, using the current pixel input value, a set of proportions and a 
set of corresponding pixel position offsets, from a family of sets of proportions and 
corresponding pixel position offsets, said set of proportions being in accordance with a 
next scanline error impulse response that approximates a function which spreads with 

5 self-convolution in proportion to a degree of self-convolution; 

(d) adding the selected proportions of the error at the current pixel to the 
neighbourhood error values of yet to be processed pixels at pixel positions offset from the 
current pixel by the selected corresponding pixel position offsets; 

and where each set of the said family of sets of proportions and corresponding 
10 pixel offsets, only includes pixel offsets corresponding to pixels on the same scanline as 
the current pixel or to pixels on the next scanline. 

According to another aspect of the invention, there is provided an apparatus for 
halftoning an image, said apparatus comprising: 

means for determining an output value of a current pixel on a current scanline 
15 using a sum of an input value for the current pixel and a neighbourhood error value at the 
current pixel; 

means for detenruning an error at the current pixel as the difference between (i) 
the sum of the input value for the current pixel and the neighbourhood error value at the 
current pixel, and (ii) the output value of the current pixel; and 
20 means for adding a proportion of the error at the current pixel to neighbourhood 

error values at as yet unprocessed pixels of a subsequent scanline in accordance with a 
next scanline error impulse response; wherein said next scanline error impulse response: 

approximates a function which spreads with self-convolution in proportion to a 
degree of self-convolution. 
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According to another aspect of the invention, there is provided an apparatus for 
halftoning an image, said apparatus comprising: 
a memory for storing a program; 

a processor for executing the program, said program comprising: 
5 code for a determining step for determining an output value of a current pixel on 

a current scanline using a sum of an input value for the current pixel and a neighbourhood 
error value at the current pixel; 

code for a determining step for determining an error at the current pixel as the 
difference between (i) the sum of the input value for the current pixel and the 
10 neighbourhood error value at the current pixel, and (ii) the output value of the current 
pixel; and 

code for an adding step for adding a proportion of the error at the current pixel to 
neighbourhood error values at as yet unprocessed pixels of a subsequent scanline in 
accordance with a next scanline error impulse response; wherein said next scanline error 
15 impulse response: 

approximates a function which spreads with self-convolution in proportion to a 
degree of self-convolution. 

According to another aspect of the invention, there is provided a computer 
program product including a computer readable medium having recorded thereon a 
20 computer program for directing a processor to execute a method for halftoning an image, 
said program comprising: 

code for a determining step for determining an output value of a current pixel on 
a current scanline using a sum of an input value for the current pixel and a neighbourhood 
error value at the current pixel; 
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code for a determining step for determining an error at the current pixel as the 
difference between (i) the sum of the input value for the current pixel and the 
neighbourhood error value at the current pixel, and (ii) the output value of the current 
pixel; and 

5 code for an adding step for adding a proportion of the error at the current pixel to 

neighbourhood error values at as yet unprocessed pixels of a subsequent scanline in 
accordance with a next scanline error impulse response; wherein said next scanline error 
impulse response: 

approximates a function which spreads with self-con volution in proportion to a 
W 10 degree of self-convolution. 

According to another aspect of the invention, there is provided a computer 
% J program for directing a processor to execute a method for halftoning an image, said 

H- program comprising: 

code for a determining step for determining an output value of a current pixel on 
O 15 a current scanline using a sum of an input value for the current pixel and a neighbourhood 

error value at the current pixel; 

code for a determining step for determining an error at the current pixel as the 
difference between (i) the sum of the input value for the current pixel and the 
neighbourhood error value at the current pixel, and (ii) the output value of the current 
20 pixel; and 

code for an adding step for adding a proportion of the error at the current pixel to 
neighbourhood error values at as yet unprocessed pixels of a subsequent scanline in 
accordance with a next scanline error impulse response; wherein said next scanline error 
impulse response: 
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approximates a function which spreads with self-convolution in proportion to a 
degree of self-convolution. 

Brief Description of the Drawings 

One or more embodiments of the present invention will now be described with 
reference to the drawings, in which: 

Fig. 1 shows an error diffusion mask as described by Floyd and Steinberg; 

Fig. 5 shows a halftone output image generated in accordance with the Floyd 
Steinberg mask shown in Fig.l ; 

Fig. 2 shows an error diffusion mask in accordance with US Patent No. 
5,353,127 (Shiau and Fan); 

Fig. 6 shows an exemplary halftone output image in accordance with bi-level 
error diffusion, using the mask shown in Fig. 2; 

Fig. 3 shows the mask positions of an exemplary error diffusion mask, designed 
in accordance with an arrangement described herein; 

Fig. 4 shows a tabular representation of the mask weight values of the exemplary 
error diffusion mask shown in Fig. 3; 

Fig. 7 shows a halftone output image produced using the mask shown in Fig. 3; 

Fig. 8 shows distribution of pixel errors for error diffusion using a single error 
line store; 

Fig. 9 shows error diffusion processing for an arbitrary pixel; 
Fig. 10 shows an alternate arrangement of error diffusion processing on a per 
pixel basis; 

Fig. 1 1 shows error diffusion processing per pixel in accordance with current and 
next scanline error impulse response functions; 

Fig. 12 provides a representation of a conservative vector field; 
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Fig, 13 shows an nth next scanline error impulse response function using an 
approximate Cauchy distribution; 

Fig. 14 shows self convolutions of a next scanline error impulse response 
function for Floyd Steinberg error diffusion; 

Fig. 15 shows self convolutions of a next scanline error impulse response 
function for US Patent No. 5,353,127 error diffusion; 

Fig. 16 shows self convolutions of a next scanline error impulse response 
function in accordance with a Cauchy distribution; 

Fig. 17 is a block diagram representation of a system for performing pixel 
processing in accordance with a first arrangement described herein; 

Fig. 19 shows an exemplary family of error diffusion masks for a second 
arrangement using a Cauchy distribution; 

Fig. 20 shows a tabular representation of a mapping between grey levels and the 
masks shown in Figure 19; 

Fig. 18 shows a block diagram representation of a system for pixel processing 
using a second arrangement of the Cauchy distribution; 

Fig. 21 is a schematic block diagram of a general purpose computer upon which 
arrangements described herein can be practiced; and 

Fig. 22 shows a block diagram of the processing per-pixel for error diffusion; 

Fig. 23 shows a block diagram of the processing per-pixel for error diffusion in 
terms of impulse response functions; 

Fig. 24 shows a 1 -dimensional view of error diffusion; 

Fig. 25 shows weights used in the Floyd-Steinberg Error Diffusion Method; 
Fig. 26 shows an influence function for Floyd-Steinberg Error Diffusion; 
Fig. 27 shows weights used in error diffusion according to Shiau and Fan; 
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Fig. 28 shows error diffusion outputs for a region of grey-level 253, for the three 
error diffusion methods indicated; 

Fig. 29 shows a wide error diffusion mask, for which the next scanline error 
impulse response approximates a Cauchy distribution; 
5 Fig. 30 shows mask weights for the various mask positions in the mask of Fig. 

29; 

Figs. 31 to 33 show next scanline error impulse response functions for Floyd- 
Steinberg, Shiau and Fan, and Cauchy methods; 

Fig. 34 shows the Cauchy distribution in terms of a vector field; 
10 Fig. 35 shows error diffusion masks for a range of widths; and 

Fig. 36 shows mappings between grey levels and masks of Fig. 35. 

Detailed Description including Best Mode 
Where reference is made in any one or more of the accompanying drawings to 
steps and/or features, which have the same reference numerals, those steps and/or features 
15 have for the purposes of this description the same function(s) or operation(s), unless the 
contrary intention appears. Further descriptive material is to be found in Appendix A. 

Some portions of the description which follows are explicitly or implicitly 
presented in terms of algorithms and symbolic representations of operations on data 
within a computer memory. These algorithmic descriptions and representations are the 
20 means used by those skilled in the data processing arts to most effectively convey the 
substance of their work to others skilled in the art. An algorithm is here, and generally, 
conceived to be a self-consistent sequence of steps leading to a desired result. The steps 
are those requiring physical manipulations of physical quantities. Usually, though not 
necessarily, these quantities take the form of electrical or magnetic signals capable of 
25 being stored, transferred, combined, compared, and otherwise manipulated. It has proven 
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convenient at times, principally for reasons of common usage, to refer to these signals as 
bits, values, elements, symbols, characters, terms, numbers, or the like. 

It should be borne in mind, however, that the above and similar terms are to be 
associated with the appropriate physical quantities and are merely convenient labels 
5 applied to these quantities. Unless specifically stated otherwise, and as apparent from the 
following, it will be appreciated that throughout the present specification, discussions 
utilizing terms such as "scanning", "calculating", "determining" "replacing", 
"generating" "initializing", "outputting", or the like, refer to the action and processes of a 
computer system, or similar electronic device, that manipulates and transforms data 
10 represented as physical (electronic) quantities within the registers and memories of the 
computer system into other data similarly represented as physical quantities within the 
computer system memories or registers or other such information storage, transmission or 
display devices. 

The present specification also discloses apparatus for performing the operations 
15 of the methods. Such apparatus may be specially constructed for the required purposes, 

or may comprise a general purpose computer or other device selectively activated or 

reconfigured by a computer program stored in the computer. The algorithms and displays 

presented herein are not inherently related to any particular computer or other apparatus. 

Various general purpose machines may be used with programs in accordance with the 
20 teachings herein. Alternatively, the construction of more specialized apparatus to perform 

the required method steps may be appropriate. The structure of a conventional general 

purpose computer will appear from the description below. 

In addition, the present specification also discloses a computer readable medium 

comprising a computer program for performing the operations of the methods. The 
25 computer readable medium is taken herein to include any transmission medium for 
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communicating the computer program between a source and a designation. The 
transmission medium may include storage devices such as magnetic or optical disks, 
memory chips, or other storage devices suitable for interfacing with a general purpose 
computer. The transmission medium may also include a hard-wired medium such as 
5 exemplified in the Internet system, or wireless medium such as exemplified in the GSM 
mobile telephone system. The computer program is not intended to be limited to any 
particular programming language and implementation thereof. It will be appreciated that 
a variety of programming languages and coding thereof may be used to implement the 
teachings of the disclosure contained herein. 
10 Modifications to error diffusion are now described, which remove worm artifacts 

in highlight and shadow regions, by design of the error diffusion mask. 

Figs. 3 and 4 show an example of an error diffusion mask, designed as described 
in this document, which generates bi-level halftone output for 8 bit per pixel source 
images which is substantially free of worm artifacts; that is, all sparse halftone patterns 
15 are well-spread. Fig. 3 shows an error diffusion mask 300 in which error for a current 
pixel 306 is distributed to neighbouring pixels on a current scanline 314, as well as pixels 
on a "next" scanline 316. Previously processed pixels are shown in full shade above a 
bold line 318, and it is noted that a current scanline 314 and a next scanline 316 contain 
pixels to which the error for the current pixel 306 is distributed. Current scanline pixels 
20 which receive error distribution from the current pixel 306 are represented by a group 
312, and next line pixels are correspondingly represented by a group 310. Fig. 4 shows a 
tabular representation of the error distribution weight values for the error diffusion mask, 
and it is noted that the table 400 is made up of four columns 402 to 408. Columns 402 
and 406 represent mask positions for the next and current scanlines respectively, and 
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columns 404 and 408 represent the proportions in which the error from the current pixel 
are distributed. 

The first arrangement is described for the case of halftoning a monochrome 8 bit 
per pixel input image to a bi-level output image. 
5 A monochrome input image with pixel values g( j being integers in the range 0 to 

255 is used to generate a bi-level halftone output image with pixel values being the 
integers 0 and 255. 

The halftone output value for each pixel is obtained by error diffusion, using an 
error diffusion mask which has been specially designed so that the next scanline error 
10 impulse response function, corresponding to the mask, approximates a Cauchy 
distribution which has been sampled and normalised. 

An example of a suitable error diffusion mask for the first arrangement is the 12- 
13 mask described in Figs. 3 and 4. This mask has weights for 12 pixel positions on the 
current scanline and 13 pixel positions on the next scanline and was prepared as described 
15 above. Use of such a specially designed error diffusion mask generates halftone output 
for which all sparse halftone patterns are well spread and substantially worm-free. 

Fig 7 shows halftone output using this mask for a source image of constant grey 
level 253. Fig. 7 shows an image 700 which is substantially free from worm artifacts, as 
exemplified by dots 702 and 704, which are well spaced from their neighbours. 
20 This description includes: 

a definition of the term "next scanline error impulse response" function; 
arguments supporting the proposition that in order for error diffusion 
with mask positions on the current and next scanlines to be worm-free, it is 
desirable that the next scanline error impulse response function should spread in 



19J201 



38l700US.doc 



-23 - 

proportion to the degree of self-convolution, or stated alternatively, should 
approximate a Cauchy distribution; 

a method for generating error diffusion masks for which the next 
scanline error impulse response function is optimised to approximate a function 
5 which spreads in proportion to the degree of self-convolution, ie. is optimised to 

approximate a Cauchy distribution. 

It will also be shown that the requirement that the next scanline eiror impulse 
response function should spread in proportion to the degree of self-convolution is closely 
related to the requirement that it should approximate a Cauchy distribution. 
10 Influence of a pixel decision on subsequent pixel decisions 

With standard error diffusion, if the modified input value of a pixel is lower than 
the threshold, then the halftone result is set low and the error distributed to neighbouring 
pixels acts to increase the likelihood that neighbouring pixels are set high. Conversely if 
the modified input value of a pixel is higher than the threshold, then the halftone result is 
15 set high and the error distributed to neighbouring pixels acts to increase the likelihood 
that neighbouring pixels are set low. 

The proportion of the (non-zero) error of an arbitrary reference pixel which 
contributes to the neighbourhood error of a subsequently processed pixel of interest, is a 
measure of a disincentive to assign the pixel of interest the same halftone output result as 
20 the reference pixel. 

Due to the typically sequential processing of pixels, and because the processing 
of each pixel includes the distribution (according to the error diffusion mask) of that 
pixel's error, being the sum of that pixel's neighbourhood error (the error distributed to 
that pixel) and that pixel's pixel-only error, a pixel contributes a proportion of its error to 
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the neighbourhood error of pixels which are beyond those to which it directly distributes 
error. 

Error diffusion imp lementable with a single error line store 
Error diffusion is now described, where: 
5 (a) pixels are processed one at a time across a scanline or row, either left to 

right or right to left, and scanlines are processed one after another from the top to 
the bottom of the image; and 

(b) the error diffusion mask includes only pixel positions in the current and 
succeeding scanline. 

10 Fig. 8 is a diagram indicating, for this class of error diffusion, the pixels to which 

error is distributed. Fig. 8 shows a representation 800 in which a current pixel 802 is 
shown on a current scanline 812, where scanline processing is from left to right as 
depicted by an arrow 806, and where previously processed pixels eg. 804 are full shaded. 
Current scanline mask positions eg. 808 denote mask values for pixel positions on the 

15 current scanline, and are designated maskc UIT [i](i>0). Next scanline mask positions eg. 
810 denote mask values for pixel positions on the next scanline, and are designated 
mask ne xi[i]. Mask positions on the current and next scanlines are identified by a 
horizontal pixel offset "i". Horizontal pixel offsets are shown in Fig. 8 by circled 
numbers 816. 

20 Error diffusion processing per pixel 

Fig. 9 shows error diffusion processing per pixel, for an arbitrary pixel of 
column i and scanline j. Fig, 9 shows a block diagram representation 900 where an input 
image grey level at pixel (ij) ie. 902 is provided to an addition process 904 which 
outputs, on a line 906, the sum of the input image grey level (gj), a total error distributed 
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directly to pixel (i,j) from pixels of the current scanline (e_curr_linei j), and a total error 
distributed directly to the pixel (ij) from pixels of a previous scanline (e_prev_lineg). 

The output signal on the line 906 is provided, on a line segment 908, to a 
threshold process 912 where it is compared against a threshold to obtain a halftone result 
5 for the pixel on a line 914, and subsequently on a line segment 916, the halftone result 
being designated ry. The modified input value on the line 906 is also provided on a line 
segment 91 0 to an adder process 920 which receives a negative value of the halftone 
result for the pixel on a line segment 918. An output from the addition process 920 is an 
error for distribution (designated e_combined ij) on a line 922. This error for distribution 

10 is provided on a line 924 to a current scanline error distribution unit (designated maskc UIT ) 
934 which distributes the combined pixel error e_combined ijs as shown by an exemplary 
arrow 932, to unprocessed pixels of the current scanline in a current line error buffer 930. 
An output from the current line error buffer 930 is provided on a line 928 (the signal 
being designated as e_curr_liney) and is delivered to an adder process 948. The 

15 combined pixel error e_combinedj d - is also provided on a line 926 to a next scanline error 
distribution unit (designated as maslwi) 936, which outputs, as exemplified by an arrow 
938, corresponding values to pixels of the next scanline which are stored in a line store 
error buffer 940. The buffer 940 provides an output on a line 942 to an error line store 
944. On a line 946, a value is retrieved from the error line store, the value being the total 

20 error distributed directly to a pixel (ij) from pixels of the previous scanline, ie. 
e_prev_linejj. 

Operation of the process of Fig. 9 is described by the following 4 steps. 
Step 1 . Derive a modified input value for the pixel as: 

gij + e_curr_lineij + ejprevjinejj (30) 
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where: gij is the input image grey level at pixel (ij), e_curr_lineij is the total error 
distributed directly to pixel (ij) from pixels of the current scanline, and e__prev_linejj is 
the total error distributed directly to pixel (i j) from pixels of the previous scanline; 

Step 2. Compare the modified input value against a threshold to obtain the 
5 halftone result for the pixel, denoted as r,j 

Step 3. Calculate an error for distribution as: 

e_combinedjj =* ey + e_currjineij + e_prev__linev,- (1) 
where: ey is the pixel-only error equal to g fj - r^ 

Step 4, Distribute the combined pixel error, e_combinedjj according to the error 
10 diffusion mask to unprocessed pixels of the current scanline and pixels of the next 
scanline. 

In step 1, the total error distributed to a pixel (ij) from previously processed 
pixels is referred to as the "neighbourhood error" at pixel (i j). The neighbourhood error, 
enbrg, modifies the pixel input value prior to thresholding and is given by: 
15 e_nbrjj = e_curr_lineij + e_prev_lineij (2) 

In step 4 error values which are portions of the combined pixel error, 
e_combinedjj are used to update error sum values in a current line error buffer and in a 
line store error buffer. 

Each error sum value in the current line error buffer is associated with a pixel 
20 position on the current scanline ahead of the current pixel - a "future pixel" of the current 
scanline; that error sum value is the sum of error values distributed directly to that future 
pixel from processed pixels of the current scanline. 

Similarly, each error sum value in the line store error buffer is associated with a 
pixel position on the next scanline - a "next scanline pixel"; that error sum value is the 
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sum of error values distributed directly to that next scanline pixel from processed pixels 
of the current scanline. 

The number of error sum values in the current line error buffer and in the line 
store error buffer need only be as large as the number of error diffusion mask positions on 
5 the current scanline and next scanline respectively. 

In step 4, once the error sum values in the line store error buffer have been 
updated, one value in the line store error buffer is complete, in that it will receive no 
further contributions, and it is transferred to the error line store. 

The number of error sum values in the error line store needs to be (and need only 
10 be) as large as the number of pixels on a scanline - hence the name "line store". 

The error sum values in the error line store generated by processing one scanline 
are used as input in the processing of the next scanline. When each pixel is processed: 

one error line store value is read (and the memory location from which 
it is read becomes available to store another error line store value) 
15 - error sum values in the line store error buffer are updated, and 

one error line store value is written. 
Fig. 10 shows an alternative implementation 1000 of the error diffusion 
processing per pixel. The input image grey level at pixel (ij), ie is provided to an 
adder process 1004. on a line 1002. The adder process 1004 outputs, on a line 1006, the 
20 sum of the input image grey level (gi d ), the total error distributed directly to pixel (ij) 
from pixels of the current scanline (ie. e_curr_liney), and the total error distributed 
directly to the pixel (ij) from pixels of the previous scanline (ie. ejprey_line;j). This 
value on the line 1006 is provided, by a line segment 1008, to a threshold process 1012 
for comparison against a threshold, which consequently produces, on a line 1014, the 
25 output result for the pixel (ie. r,.j) on a line 1016. 
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An inverse of the halftone result (ie. negative r;.j) is provided on a line 101 8 to an 
addition process 1020, along with the modified input value 
(ie. gjj + e_curr_linei d + e_prev_linei.j) on a line 1010. An output from the addition 
process 1020 is provided on a line 1022, this being the error for distribution 
5 e^combined^. This error for distribution is provided on a line 1024 to the current 
scanline error distribution unit (designated as mask^) 1028, which outputs, as 
exemplified by an arrow 1030, error distributions to unprocessed pixels of the current 
scanline in a current line error buffer 1032. The buffer 1032 consequently outputs the 
total error distributed directly to pixel (ij) from pixels of the current scanline (ie 
nj 10 e_currjine,j) on a line 1034 to an addition process 1036. The error for distribution on the 

m line 1022 is also distributed, on a line 1026, to an error line store 1038. From the error 

SJ line store a weighted sum of error values is retrieved on a line 1040 to a next scanline 

\=-% error distribution unit (designated as mask„ cxt ) 1042, which distributes the weighted sum, 

Hj as exemplified by an arrow 1044, to a line store error buffer 1046. The buffer 1046 

□ 15 provides, on a line 1048, the total error distributed directly to the pixel (ij) from pixels of 

the previous scanline (ie. e_prev_linejj) to the addition process 1036. 

The difference between the 2 implementations, is that in the alternative imple- 
mentation: 

(a) the combined error at a pixel is written to a location in the error line 
20 store; 

(b) the contribution to neighbourhood error by direct distribution from 
pixels of a previous scanline is determined as a weighted sum of error values 
retrieved from the error line store. That is, in the alternative implementation, 
distribution of error to pixels of the next scanline is achieved by gathering error 
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values from the line store, rather than distributing error amongst line store 
locations. 

An analysis is presented below concerning self-convolution of a next scanline 
error impulse response function. That analysis is described referring to Fig. 9; however 
the same analysis also applies to the alternative implementation of Fig. 10. 
Next scan line error impulse response function for error diffusion im plementable with a 
single erro r line store 

Error diffusion is now considered for the case where 100% of the combined pixel 
error is distributed to as yet unprocessed pixels (step 4 above), and the fraction of error 
distributed to the next scanline is considered to be positive and non-zero. 

Further, error diffusion is considered for the case where the error diffusion mask 
weights do not change from pixel to pixel. That is, it is assumed that the processing per 
pixel is spatially invariant. 

The distribution of the combined pixel error of an arbitrary reference pixel (but 
which is not near the left and right edges of the image) is now considered according to the 
error diffusion mask, Part of that combined pixel error is distributed to pixels of the next 
scanline, and the remainder of the combined pixel error is distributed to pixels of the 
current scanline. This remainder error remains subject to error distribution according to 
the error diffusion mask as part of the processing of the current scanline. With the 
complete processing of the current scanline, effectively all of the combined error of the 
reference pixel is distributed to pixels of the next scanline. 

The distribution of combined pixel error of an arbitrary reference pixel of the 
current scanline to pixels of the next scanline, as a result of the complete processing of the 
scanline, defines a "next scanline error impulse response" function which is denoted as 
hnexi- To be precise, h nwc is defined as a function which maps integers to real values; 
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where h ne xt[i] is that fraction of the error at a reference pixel which, following the 
complete processing of the scanline of the reference pixel, is distributed to the pixel on 
the next scanline horizontally offset i from the reference pixel. This is shown 
illustratively with reference to Fig. 13. 

Self-convolutions of the next scanline error impulse response function 

As the domain of the function hnexiis the set of integers, it can also be considered 

as a (two-sided) sequence. 

The set of input image grey levels, g; J5 is considered for a scanline j as a 

sequence, and the sequence is denoted as g. 
Similarly: 

the sequence of error sum values, e_prev_linejj, being sums of errors 
distributed directly to pixels of scanline j from scanline j-1, is denoted as 
e_prev_linej; 

the sequence of halftone output image values for scanline j is denoted as 

the sequence of pixel-only errors for scanline j is denoted as ej. 
For the class of error diffusion considered, the processing performed for each 
scanline, j, can be considered to take as input the following: 

ejpreyjinej, gj (31) 
and produces as output the following: 

r,-, e^ e_prev_Hnej+i (32) 
For the implementation of Fig. 9, as part of the processing of scanline j, the 
sequence of values e_prev_linej are read from the error line store and the sequence of 
values e_prev_Iine,+i are written to the error line store. 
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From the definition of the next scanline error impulse response function, the 
following can be written: 

e_prevjine j+ i = (ej + e_prev_linej) * Iw, (3) 
That is, the sequence of sums of error values distributed directly to pixels of 
5 scanline j+1 from scanline j can be represented as the sum of the error sequence for 
scanline j and the sequence of sums of error values distributed directly to pixels of 
scanline j from scanline j-1 convolved with the next scanline error impulse response 
function. 

Here the convolution of 2 sequences is denoted by and the sequence formed 
10 by the convolution of 2 sequences, f and g is defined as: 

(f*g)[i] = Z keZ f[i-k]g[k] (33) 
where Z is the set of all integers. 

Assuming that the next scanline error impulse response function is the same for 
each scanline, (3) can be applied recursively, and the sequence of sums of error values 
15 distributed directly to pixels of a scanline from the previous scanline can be written as a 
weighted sum of pixel only errors of pixels of preceding scanlines. 
That is, 

e_prev_linej => (e^ + e_preyjine>0 * hncxt (33)* 
= (ej-i + (q-2 + e _prev_linej. 2 ) * h nex t) * Iw (34) 

20 and so on, giving 

e_prev_linej = 2, <j e, * h next * (J - ,) (4) 
where hnext*" denotes the sequence formed as the convolution of the sequence h ne xc with 
itself 'n- 1 ' times. So that 

hnext*' =h ncxl (35) 
25 hnext* 2 " hnext * ^nem and so on. (36) 
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The current scanline er r or imp ulse response function in terms of the error diffusion mask 

The error diffusion mask can be represented using 2 functions which map the 
integers into real values, mask^rr and mask^t, where: 

maskcurrfi] is the error diffusion mask weight for that pixel on the current 
scanline horizontally offset by T from the reference pixel being processed; and 

mask„ CK i[i] is the error diffusion mask weight for that pixel on the next scanline 
horizontally offset by T from the reference pixel being processed. 

Fig. 8 shows horizontal pixel offsets 816 from a reference pixel, together with 
mask positions 808 and 810 on the current and next scanline respectively. 

Due to the sequential processing of pixels of the current scanline, 

mask^i] = 0 for i <= 0. (37) 

As a consequence of the error diffusion processing per pixel and from 
considering fig. 9 it is seen that: 

e_cuTT_linei (i ~Z keZ e_combined; . kj . mask cun -[k] (38) 

That is, the corresponding sequences are related by the convolution operation: 

e_curr_linej = e_combinedj * maskc UTT (5) 

The sequence e_combinedj is the same as (ej + e_currjinej + e_prev_linej), and 

so: 

e_cmr_linej = (ej + e_curr_linej + e_prev_linej) * mask^ (6) 
By repeatedly replacing e_curr_linef in the right hand side of (6) with the entire 
right hand side of (6), the following relationship results: 

e_cuiT_linej = (ej + e_prev_linej) * (masko urT + maskc urT *maskcun- + ••) (39) 
That is, 

e_curr_linej = (ej + e_prev_linej) * h cun (7) 
where h ewrr is the function (sequence) given by: 
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hcun- = maskcurr + maskcun-* 2 + maslw* 3 + . . (8) 
hcun- is equivalently defined by the recursive definition: 

hcurr = (8 + hcurr) * maskc urr (9) 
where 6 is the delta function (sequence) given by: 

8[k] = lfork = 0 (40) 
5[k] - 0 for k * 0 (41) 
hcun- is referred to as the current scanline error impulse response function. 
The next scanline error impulse response function in terms of the error diffusion mask 

As is the case for the current scanline error impulse response function, the next 
scanline error impulse response function is determined by the error diffusion mask. 

As a consequence of the error diffusion processing per pixel and from 
considering Fig. 9 the following relationship can be expressed: 

e.jireyjinei, j+1 = S k s z ecombined; . Kj . maskn CXl [k] (42) 
That is, the corresponding sequences are related by the convolution operation: 

e_jprev_linej+i = e_combinedj * maskn ext (10) 
As above, using the fact that the sequence e_combinedj is the same as 
(ej + e_currjine, + e_prev_linej), the following is seen: 

e_prey_lme-j+i (ej + e_curr_linej + e_prev_linej) * mask„ cw (1 1) 

By replacing e_curr_hnej in the right hand side of (1 1) with the right hand side of 
(7), the following is seen: 

ejH-eyJinej+i - (e,- + e_prev_linej + (ej + e_prev_linej) * Iw) * masknext (43) 
which can be re-written as 

e_prev_line j+ i = (ej + ejpreyjinej) * ((5 + Iw) * mask^xO (44) 
Comparing this with (3), it is seen that that the next scanline error impulse 
response function is given by: 
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hnext - (5 + hc ur r) * mask ncx; (12) 
By substituting the expansion of iw of (8) into (12) it is seen that 
rw - masknext * (5 + masker + mask^* 2 + maskc Urr * 3 + . . ) (13) 
Current and next scanline error impulse response functions for standard error diffusion 

To assist an understanding of the above definitions of the current and next 
scanline error impulse response functions, the current and next scanline error impulse 
response for standard (Floyd Steinberg) error diffusion are now described. 

For Floyd Steinberg error diffusion, the error diffusion mask is defined by the 
following mask functions: 

mask culT [i] = Ofori*l 1 

mask 01I1I [l] = 7/16 J (45) 



mask nexI [i] = Ofbri<-landi>l ] 
mask nOT [-l] = 3/16 ! 
mask ncxt [0J = 5/16 \ 
mask ^[1] = 1/1 6 J 



(46) 



For Floyd Steinberg error diffusion, the current scanline error impulse response 
function is given by: 

h CUIT [i] = 0fori 



ri<=0 1 
6/ fori>=lJ 



(47) 

h tlflT [i] = (7/16)' - ' 
which is consistent with (8) and (9). 

For Floyd Steinberg error diffusion, the next scanline error impulse response 
function is given by: 



h Mxl [i] = 0fori<-l 1 
h llcxt [-l] = 3/16 | 
h»«i[0] = 5/1 6 + 7/1 6* 3/16 

h next [i] = (l/16 + 7/16*5/16 + (7/16) 2 *3/16)*(7/16) il fori>0 J 



(48) 
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which is consistent with (12). 

Description of error diffusion processing in terms of the current and next 
scanline error impulse response functions 

This description provides, firstly, further clarification of the current and next 
5 scanline error impulse response functions, and secondly, it is the basis for later discussion 
concerning desirable characteristics of the current and next scanline error impulse 
response functions. 

Fig. 1 1 shows error diffusion processing per pixel in terms of the current and next 
scanline error impulse response functions, hc uri - and Kan- The input image grey level at 

10 pixel (ij), ie. g^, is provided, as depicted by an arrow 1102, to an addition process 1104. 
The addition process 1 104 provides, as depicted by a line segment 1106, the sum of the 
input image grey level (gij), and the total error distributed directly to the pixel (ij) from 
pixels of the previous scanline (ie. e__prev_lineij) on a line 1108, to an addition process 
1 1 12, and also, on a line 1 1 10, to an addition process 1 124. The addition process 1112 

15 provides an output on a line 1 114 to a threshold process 1116, which compares the value 
against a threshold, consequently outputting the halftone result for the pixel on a line 
1118, An inverted value of the halftone result (ie. -r id ) is provided on a line 1 122 to the 
addition process 1124. 

The addition process 1124 outputs, on a line 1126, and on a line segment 1128, an 

20 aggregate value to a current scanline error impulse response distribution unit 1134 
(designated h* wr ) which provides, as exemplified by an arrow 1134, values distributed to 
pixels of the current scanline in a buffer 1 136. From a buffer 1 136 on a line 1 138, a value 
labelled as e_curr_linejj is provided to the addition process 1112. In a similar manner, the 
signal on the line 1126 is provided, on a line segment 1130, to a next scanline error 

25 impulse response distribution unit (designated h next ) 1 140 which, as exemplified by an 
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arrow 1142, distributes a value in accordance with the next scanline error impulse 
response h next to pixels of the next scanline in a line store error buffer 1 144. The buffer 
1 144 consequently outputs, on a line 1 146, a value to the error line store 1 148. From the 
error line store, a value labelled as e_prev_Hnejj is retrieved and provided to the addition 
process 1104. 

Now the labels, e_curr Jinejj and e.jsrevjiney in Fig. 1 1 are consistent with the 
previous definitions of these terms. This is because 

1 . the sequence resulting from applying the current scanline error impulse 
response to the sequence ej + e_prey_linej is the same as the sequence resulting 
from applying the current scanline mask coefficients to the sequence 
e_combinedj, which can be seen from (5) and (7) follows: 

e_curr_linej = e_combinedj * masker = (ej + e_prev_linej) * hcurr (49) 

2. the sequence resulting from applying the next scanline error impulse 
response to the sequence ej + e j>rev_linej is the same as the sequence resulting 
from applying the next scanline mask coefficients to the sequence e_combinedj, 
which can be seen from (3) and (10) as follows: 

e_prey_line j+ i = e_combinedj * raaskn CXt = (ej + e_prev_linej) * h ne xt (50) 
That is, if hc UTT and hn ext are defined in terms of maskc UIT and mask nex t according to 
(8) or (9) and (12) or (13), then the processing described in relation to Fig. 11 is 
equivalent to the processing described in relation to Figs. 9 and 10 in that they produce 
the same output. In each case, for each pixel, the same quantity, gij + e_curr_line;,j + 
e__prev_line,j is applied to the threshold operation. 

In general, for implementations of error diffusion, the processing models of Figs. 
9 or 10 are preferred over the processing model of Fig. 11. This is because error diffusion 
masks with non-zero coefficients at a small number of pixel offset positions, correspond 



191201 



581700US.doc 



-37- 

to current and next scanline error impulse response functions which typically have an 
unbounded number of non-zero function values. 

However, the current and next scanline error impulse response functions and the 
processing model of Fig. 1 1 are useful for describing desirable features of error diffusion 
5 processing, and also for performing experiments to establish desirable features of error 
diffusion processing. 

The neighbourhood error at a current pixel, is the total error distributed directly 
from previously processed pixels to the current pixel, and is used to modify the current 
pixel's input grey value prior to thresholding. Considering the sequence of 
10 neighbourhood errors for a scanline, the following is seen: 

e_nbrj = e^currjinej + e_prev_linej (50)* 

= (e, - + e_prevjinej) * rw + (ej-i + e_prev_linej.i) * rwt (5 1) 

giving, 

enbtj = Ej * lw + Eh * h next (14) 
15 The above equation (14) neatly isolates the contribution to neighbourhood error 

from pixels of the current scanline and pixels of the previous scanline in terms of a 
common error value per pixel, That common error value per pixel warrants a special 
name. The (previous scanlines) "modified pixel error" for a pixel (ij) Ej j is defined as: 

Ejj = ey + e_prev_linejj (52) 
20 Equation 14 provides assistance in inferring some conclusions regarding 

desirable characteristics of error diffusion. 

Desirable attributes of thf - ri^vt scanl ine error impulse res ponse function: requirementof 
unity sum 



191201 



5817O0US.doc 



-38- 

As noted earlier, in order that the average halftone image output matches the 
average image input, it is a requirement that the sum of the error diffusion coefficients 
should be 1. That is, 

Sj maskcun.fi] + Zj mask ncM [i] = 1 (25) 
It follows from equations 13 and 15, assuming that the absolute value of the sum 
of error diffusion mask coefficients on the current scanline is less than 1, that the sum of 
coefficients of the next scanline error impulse response function is 1 : 

2ilw[i] = l (16) 
In summary, for the processing model of Figs. 9 and 10, involving the 2 
functions mask^n- and mask,,^; the requirement that average halftone image output 
matches average image input is summarised by (15), whereas for the processing model of 
Fig. 11, the same requirement is summarised by equation 16, involving only one function 

hnext- 

That is, in the processing model of Fig. 11, the current scanline error impulse 
response function can be varied independently, without compromising the requirement for 
average halftone output to match average image input. 

Desirable attributes of the next s canline error impulse response function: requirement of 
left - right symmetry 

The next scanline error impulse response function as defined by equation 3, and 
as depicted in Fig. 11, does not show any dependence on whether current scanline 
processing is left to right or right to left. 

If it is assumed that current scanline processing does a very good job of setting 
halftone output values for the scanline, and consequently generating pixel only error 
values for the scanline, then it is reasonable to assume that subsequent processing of the 
next scanline could be left to right or right to left and does not need to compensate for any 
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left-right bias in the pixel only errors of the current scanline. It seems clear then that the 

next scanline error impulse response function should not be biased to left or right but be 

left-right symmetric. 

A requirement is thus defined that the next scanline error impulse response 
5 function should be left-right symmetric. 

In fact, some signal shifting or phase distortion due to the current scanline error 
impulse response function is inevitable because it is a causal filter. However, the 
signal shifting and phase distortion can be made small. So, it is reasonable to 
ignore compensation for left-right bias by use of the next scanline error impulse 
10 response function, at least for discussion of appropriate spreading of the next 

scanline error impulse response function. 

Desirable attributes of the current and next scanline error impulse response functions: 

requirement of being monotonic decreasing with increasing horizontal pixel offset 

The current scanline error impulse response function measures the degree to 
15 which the halftone error at a current pixel on the current scanline should be taken account 

of in the halftone decision of a subsequently processed pixel on the current scanline. 

Also, having concluded, at least to a first approximation, that the next scanline 

error impulse response function should be left-right symmetric, there is no preference in 

whether a scanline is processed in the same direction as the preceding scanline. The next 
20 scanline error impulse response function contributes to the amount to which the halftone 

error at a current pixel on the current scanline is taken account of in the halftone decision 

of a pixel on the next scanline. 

It is further suggested that both the current and next scanline impulse response 

functions should decrease monotonically with increasing horizontal pixel separation. 
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Desirable att ributes of the next scanline error impulse response function: requirement of 
similar shape of current and next scanline error impulse response functions 

Having concluded, at least to a first approximation, that the next scanline error 
impulse response function should be left-right symmetric, it can be seen using (52) and 
5 (4) that each of the modified pixel error values, Eij, is a weighted sum of pixel-only errors 
and is obtained by convolution of pixel-only scanline sequences (ej), with left-right 
symmetric function h ncxt *(j-i)- 

Thus, (2) and (14) isolate the asymmetric and symmetric influence of pixel only 
error of previously processed pixels on the halftone decision at a current pixel. 
10 That is, the term 

e_curr_lmejj = (Ej * Iw)[i] (53) 
is an asymmetric contribution to the neighbourhood error of pixel (i j), being the 
error distributed to the pixel directly from pixels of the current scanline; whereas the term 
e_prevjineij = (Ej.! * h ncJtt )[i] (54) 
15 is a symmetric contribution, being the error distributed to the pixel directly from 

pixels of the previous scanline. 

In a similar argument to that for inferring monotonicity and left-right symmetry, 
it is suggested that the contribution to the neighbourhood error of the current pixel, from 
the modified error of previously processed pixels of the current scanline, should reduce 
20 with increasing horizontal separation in a similar fashion to the contribution from the 
modified error of previously processed pixels of the previous scanline. In this way no 
particular horizontal separation is favoured on either the current or previous scanline. 
Desirable attributes of the next scanline error impulse response function: requirement of 
appropriate spreading 
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Experiments have been performed, observing halftone output when the next 
scanline error impulse response function has various shapes. In each case, the current 
scanline response function has the same shape as the next scanline response function, 
being determined from the next scanline response function by multiplication by a step 
function. Note that in these simulations the number of non-zero values used for the 
current and next scanline response functions was set high to avoid misleading results due 
to truncation of the functions. 
Experiment 1 

In this case the next scanline error impulse response function is a sampling of the 
2-sided exponential distribution. 

h CUT [i] = a.w i fori>=l 
b ue . l [i3 = ((l-w)/(l + w)).w |i| 

where a is a positive constant and w is a constant with 0 < w < 1 

Note that with a = 1, the above current and next scanline error impulse response 

functions correspond to the following mask functions 

mask CUJT [ 1 ] = w and mask cm . [i] = 0 for i * 1 
mask^M^- w ) 2 w '" fori <=0andmask DCXl [i] = 0 fori >0 

With a = 1 and w = 1/2, the above mask functions and functions h wn and h nex; 
provide an extension to US Patent 5,353, 127 (Shiau & Fan) where the extended 
distribution set of the error diffusion mask is left extended without limit. 
Observations for experiment 1 

When a = 1 and w = 1/2, halftone output at very low or very high grey levels 
suffers from worm type artifacts where the horizontal separation of minority pixels is too 
small. These artifacts can be reduced by increasing either of the parameters a or w. 
However modifying these parameters so that minority pixels at very low or very high 
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grey levels are better separated, introduces other unpleasant artifacts, including artifacts 
where for other grey levels the vertical separation between pixels becomes too small. 
Experiment 2 

In this case the next scanline error impulse response function is a sampling of the 
gaussian distribution. 

kcur W = a • exp(~ n . (i / b) 2 > for i >= 1 
hn CXt n = cexp(-^.(i/b) 2 ) 

where a, b, c are positive constants, with c chosen so that Li h ncxl [i] = 1 
The parameter b defines the width of the gaussian distribution. 
Observations for experiment 2 

When a — 1/2 and b = 4, halftone output at low or high grey levels suffers from 
worm type artifacts where the horizontal separation of minority pixels is too small. These 
artifacts can be reduced by increasing a or b. However, again, modifying these parameters 
so that minority pixels at low or high grey levels are better separated, introduces other 
unpleasant artifacts, including artifacts where for other grey levels the vertical separation 
between pixels becomes too small. 
Experiment 3 

In this case the next scanline error impulse response function is a sampling of the 
Cauchy distribution (also known as the Lorentz distribution). 

kcurrW = a.l/(b 2 + i 2 )fori>=l 
h^W-c.l/Cb 2 +i 2 ) 

where a, b, c arc positive constants, with c chosen so that 2; h ncxt [i] = 1 

The parameter b controls the spread of the distribution. 
Observations for experiment 3 
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When a = 1 and b = 1, the halftone patterns for all grey levels with sparse 
halftone patterns are well spread. Also the general quality of the halftone output is high 
across all grey levels. 

With S bit monochrome halftoning, the extreme sparse halftone patterns 
correspond to grey levels 1 and 254. By allowing fractional grey levels less than 1 the 
behaviour of error diffusion can be observed for halftone patterns which are much more 
sparse than the halftone patterns for grey level 1. 

It is noted that the very sparse halftone patterns for sub-unity fractional grey 
values are also well spread when the Cauchy distribution is used for the current and next 
scanline response functions. That is, error diffusion processing using the Cauchy 
distribution appears capable of generating well spread halftone patterns no matter how 
sparse the patterns. 

When b is increased well above 1 the vertical separation of minority pixels in 
sparse halftone patterns is too small; when b is decreased well below 1 the horizontal 
separation of minority pixels in sparse halftone patterns is too small. 
Experiment 4 

In this case the next scanline error impulse response function is a generalisation 
of the Cauchy or Lorentz distribution. 

h eutT [i] = a .l/(b 2 +i z ) p fori>=l 
K cxt [i] = cA/(b 2 + i 2 y 

where a, b, c, p are positive constants, with c chosen so that Ej h nex t[i] = 1 
Observatio ns for experiment 4 

With p set to be greater than 1, the observations are similar in character to those 
for experiments 1 and 2. 
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With p set to be less than 1 (but necessarily greater than 1/2), the halftone 
patterns for low or very high grey values are well spread horizontally, but the horizontal 
separation for intermediate grey values is too small. 
Conclusion from the experiments 

As noted, the Cauchy distribution of experiment 3 generates well spread halftone 
patterns across the entire range of grey values with sparse patterns, while the other 
distributions fail to do so. It is believed the reason for the constrasting behaviours lies in 
the 'spread' of the distributions and in the spread of the self-convolutions of the 
distributions. 

Expressing neighbourhood error in terms of the current and next scanline error 
impulse response functions and pixel only errors, from (2),( 4), and (7) it is seen that: 

e_nbrj - e_curr_line, + e_prev_linej (59)* 
= (e,- + e_prev_Jinej) * hc Un - + e_prev_lmej (60) 
= (ej + 2, <j e, * hW 0 - 1 *) * 1W + Si <j q * hnext* 0 " 0 (61) 
Accordingly, if the strength of the next scanline error impulse response function 
and its self convolutions is weak beyond some width, then insufficient account is taken of 
pixel only error at pixels horizontally separated by that width or more. Hence the 
disincentive to place minority pixel results at that width is insufficient, leading to sparse 
patterns which are not spread enough horizontally. 
Self-convolution of the Cauchy distribution 
The Cauchy distribution is given by 

f(x)=(l/*)(b/(b 2 + x 2 )) (17) 
where x is a real number and b is a real positive constant. 

By change of variable from x to 0, where tan 8 == b / x, it can be shown that, 
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The Cauchy distribution satisfies the following self-convolution equations: 
(1/tx) (b / (b 2 + x 2 )) * (1/tt) (b / (b 2 + x 2 )) = (1/tc) (2b / ((2b) 2 + x 2 )) (63) 
and for n - 1, 1, 3 .. 

(1/tc) (b / (b 2 + x 2 )) * n = (1/tt) (nb / ((nb) 2 + x 2 )) (1 8) 

These results can be established from Fourier Transform theory. The following 
definitions and notation for the Continuous Space Fourier Transform are used: 
Analysis equation / forward transform. F(w) - f(x) e" jvot dx (64) 

Synthesis equation / reverse transform f(x) = (l/27t) J"^ F(w) e iwx dw (65) 

Fourier Transform pair: f(x) <=> F(w) (66) 

The Fourier Transform pair corresponding to the double exponential distribution, 
with 'a' being a positive real constant, can be shown to be 

e- aW <i>2a/(a 2 + w 2 ) (67) 
By the symmetry / duality property of the Fourier Transform 

( 1/tc) a / (a 2 + x 2 ) <^ e- aN (68) 
By the convolution property of the Fourier Transform, with b = a, 

((1/7!) b / (b 2 + x 2 )) * n o e nb|w| (69) 
From this, (IS) can be deduced. It can be stated that the Cauchy distribution is 
determined by its Continuous Space Fourier Transform which is a two-sided exponential 
function. 

Discrete convolution 

Related equations exist for discrete convolution, 
Fork eZ, let 

h[k] - (b / w) ((1 + (-l)"" 1 e" bn ) / (b 2 + k 2 )) (19) 

then for n - 1, 2, 3 .. 

h[k] * n = (nb / ri) (( 1 + (- 1 ) w e nb,c ) / ((nb) 2 + k 2 )) (20) 
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The above equations can be established using Fourier Transform theory, either using 
Fourier Series or using the Discrete Space Fourier Transform. The following notation is 
5 used for the Discrete Space Fourier Transform: 

Synthesis equation: h[k] =(l/27t) \ 2w H^*) e jwk dw (70) 

Analysis equation: H(e jw ) = £ k mZ h[kj e" Jwk (7 1 ) 

the use of the argument, e* w , rather than w, indicates that the function, H, is periodic in w, 
with period In, so that: 

10 H(e i(,v+2 ^ = H(e 3 ' w ) for all w. (72) 

The sequence (or discrete signal or discrete space function), h[k] of (19), has the 
Discrete Space Fourier Transform given by: 

H (e ,w ) - e" bW for -n < w ^ n . (73) 
It can be stated that the discrete space function of (19) is determined by its 
15 Discrete Space Fourier Transform which is a replicated two-sided exponential function. 

Equation (20) follows by considering the self-convolution of the function of 
equation 19 and from the convolution property for the Discrete Space Fourier Transform, 
(That is, the convolution of 2 discrete time functions is the product of their Fourier 
Transforms.) 

20 Equation (21) follows by considering the value of H at w = 0. 

The sequence, or discrete space function, 

h[k] = (sinh(b7t) b / tv) (- 1 ) k / (b 2 + k 2 ) (74) 
has the Discrete Space Fourier Transform given by 

H(0 = coshfbw) = (e bw + e' bw ) 12 for - tc <! w < tz (75) 
25 Considering the value of H at w = n gives 
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2 k6Z (b/7t)(l/(b 2 + k 2 )) = coth(b7r) (22) 
Considering equations (19) and (22) with b - 1, it is seen that e"" is much less 
than I, being approximately 0.0432 and tanh(7i) = 1 / coth(7i) is close to 1, being 
approximately 0.9963. 

So, defining h[k] according to the following equation: 

h[kj = (tanh(7c) / 7i) (1 / (1 + k 2 )) (23) 
it is deduced by setting b=l in (22) that: 

S k€Z h[k] = l (76) 
Further, h[k] according to (23) satisfies the following approximation: 

h[k]«(l /7t)(l/(l +k 2 )) (77) 

and 

h[k]* n * {(1 / tc) (1 + (-ly^O / (I + k 2 )}*" - (n / ti) (1 + (-l^-V™) / (n 2 + k 2 ) (78) 
« (n / 7t) (1 / (n 2 + k 2 )) forn-1,2,3 .. . (79) 
The function of equation (23) can be considered a discrete approximation to the 
(continuous) Cauchy distribution with 'b' parameter equal to 1; it is derived from that 
Cauchy distribution by sampling and by normalising so that the sum of the function 
values is 1. 

The Cauchy distribu tion spreads in proportion to the degree of self-convolution 

Considering (18), it is seen that the Cauchy distribution has the remarkable 
property that repeated self-convolution preserves the form of the distribution. The self 
convolution of degree n of the Cauchy distribution is a scaled copy of the original 
distribution on a scaled axis: 

f(x) * n - (1/ti) (nb / ((nb) 2 + x 2 )) = (1/n) (Viz) (b / (b 2 + (x/n) 2 )) (80) 

So that 

£(x) * n = (1/n) f(x/n) (24) 
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The variance of the Cauchy distribution is not finite. An alternative measure of 
the width of the distribution is provided by the "equivalent width" as described in 'The 
Fourier Transform and its Applications" by R.N. Bracewell (page 148) as follows. 

The equivalent width of a function is the width of the rectangle whose height is 
5 equal to the central ordinate and whose area is the same as that of the function: 

mr (81) 

From (17) and (18), or from (24) it can be seen that the central ordinate of the 
self-convolution of degree n of the Cauchy distribution is 1/n times the central ordinate, 
f(0) 3 of the Cauchy distribution. The area under f(x) is 1 ; and as a consequence the area 
10 under f(x)* n is also 1. So the equivalent width of the Cauchy distribution increases in 
proportion to the degree of self-convolution. 

The flux density of a point source, conservative vector field through a straight line is 
described by th e Cauchv distribution 

The manner in which the Cauchy distribution spreads in proportion to the degree 
15 of self-convolution can also be shown graphically. 

Fig. 12 shows a vector field, A, with vector magnitude given by 1/r where r is the 
distance from a source point, O (ie. 1216), and vector direction pointing away from O. 
This vector field is conservative in that the flux out of any simple closed curve, not 
enclosing the source point, is zero. The magnitude of the vector field can be considered 
20 as a signal which spreads uniformly in 2 dimensions out from O, preserving its strength - 
that is, reducing in proportion to the length of its circular wavefront. 

Considering the flux across a horizontal line, L (ie. 1202), which has P (ie. 1214), its 
nearest point to O (ie. 1216), at a distance b (ie. 1210) from O (ie. 1216), the flux density 
function, f(x), is defined as a function of the offset, x, along the line from P, such that 
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i 

J*(x)<3x (82) 

is the flux of the vector field through a segment [s,t] of line L. 

Let the angle which a point at offset x along the line L, makes with the line from 
O to P be 6 (ie. 1220), then fi» is given by 

f(x) = c (1/r) cos e (83) 

with c being a constant. 

Choosing c to be 1/ti so that the total flux through the line L is 1, we have 

f(x) - (l/7t) (1/r) (b/r) = (1/tt) b / (b 2 + x 2 ), (84) 
which is the Cauchy distribution. 

The offset, x, and angle, 6, are related by: 

x - b tan8 (25) 
dx/de = b (1 + tan 2 G) (85) 
d9/dx = b/(b 2 +x 2 ) (26) 
Using (25) and (26), confirms that the flux through a segment [s,t] (ie. 1218) is 
proportional to the angle, 0 SI (ie. 1224), subtended at O as shown in the following: 

Considering another line, L' (ie. 1204), with nearest point at a distance nb (ie. 
1206) from O and which is parallel to L. Segment [ns, nt] (ie. 1222) is the projection of 
segment [s,t] (ie. 1218) onto the line L" (ie. 1204), and subtends the same angle 9 st (ie. 
1224) at O (ie. 1216). 

The flux through segment [ns, nt] (ie. 1222) is also equal to 9 st /7* and is 
associated with a flux density function for line L' equal to the self-convolution of f(x) of 
degree n as follows: 
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r m^dx = f nt If(^)dx = f nt I 2* dx = f *i ffuMu - f8tt 

Jns J„ n V ) n s 7r(inb y + ^ ) ax -i s ^ tWaxX -~ ( 8? ) 

So It is seen that flux density functions given by the Cauchy distribution and its 
self convolutions correspond to a signal radiating from a point source which preserves its 
strength within any wedge formed by 2 rays emanating from the signal source point, 
Desirability of the next scanline error impulse response function snreadine in proportion 
to the degree of self-convolution and approximating a Cauchv distribution 

It has been shown that the spread of the next scanline error impulse response 
function and its self convolutions relates to the horizontal separation between minority 
pixels in sparse halftone patterns. 

It has been shown by experiment that when the next scanline error impulse 
response function samples a Cauchy distribution with the a and b parameters close to 1, 
sparse halftone patterns are well spread and worm-free. As well, when the next scanline 
error impulse response function is derived from equation 19, sparse halftone patterns are 
again well spread and worm-free. 

It has also been shown that self convolutions of the Cauchy distribution spread in 
proportion to the degree of self-convolution and maintain a radially uniform distribution 
of signal strength - neither shifting signal strength to the centre nor shifting it towards the 
extremities. 

A real, symmetric distribution which spreads in proportion to the degree of self- 
convolution must be a Cauchv distribution 

It can be shown that a real, symmetric continuous impulse response function 
which preserves the magnitude of the impulse, spreads in proportion to the degree of self- 
convolution and which acts as a smooth (continuous) low pass filter, dampening all non- 
zero frequencies must be a Cauchy distribution. 

Consider a positive real valued impulse response function f(x). 
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The condition that f(x) spreads in proportion to the degree of self-convolution 
can be written as follows: 

f(x) * n = a f(x/n) for some constant a (88) 

Let F(w) be the Fourier Transform of f(x). The condition that f(x) preserves the 
magnitude of the impulse can be written as follows: 

Fco)=r. 00 fl:x)dx = i (89) 

By reference to the Fourier convolution theorem, the condition that the 
magnitude of the impulse is preserved by self-convolution can be written as follows: 

f%, fix)*" dx = [F(0)] n = 1 (90) 
Also, the integral of f(x/n) over all real values is n. So it can be deduced that 

a=l/n (91) 
So f(x) satisfies (24) as follows: 

fix) *" - (1/n) fix/n) (92) 
Taking the Fourier Transform of each side of (24) gives the following: 

(F(w)) n = F(nw) for n = 1 , 2, 3 (93) 
Differentiating the above equation with respect to w gives the following: 

n . (F(w)) n " 1 . F'(w) = n . F'(nw) (94) 
By dividing (93) by (94), the following results: 

F(w) / F(w) =F(nw) / F'(nw) for n = 1 , 2, 3 . . (95) 
By the condition that F is continuous (smoothly dampening all non-zero 
frequencies) it can be deduced that: 

for w >= 0, F(nw) / F'(nw) = ci, a constant (96) 
Solving this differential equation gives the following: 

F(w) = c 2 e B w for w > 0 for constants c 2 and B (97) 
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Because f(x) is real-valued and symmetric, F(w) must be also; so that C2 and B 
are both real. For F(0) - 1, c 2 = 1. For f(x) to act as a low pass filter, B must be negative. 
From the above, it can be deduced that: 

F(w) = e" b|w| for a positive constant b. (98) 
5 Using the inverse Fourier Transform gives the following: 

f(x) = b / (b 2 +x 2 ) as claimed. (99) 

Cauchv error diffusion 

The (continuous) Cauchy distribution spreads in proportion to the degree of self- 
convolution. It is believed that this spreading property is desirable for the (discrete) next 

10 scanline error impulse response function. Moreover, it is envisaged that this spreading 
property is desirable for error diffusion irrespective of the number of scanlines with mask 
positions, and is desirable for other neighbourhood halftoning techniques. Similarly, it is 
believed that it is desirable that, in some sense, the (discrete) next scanline error impulse 
response function should approximate a (continuous) Cauchy distribution. 

15 The name "Cauchy error diffusion" is defined to be error diffusion where the 

next scanline error impulse response function is designed so that it spreads approximately 
in proportion to the degree of self convolution or where the next scanline error impulse 
response function is designed to approximate a Cauchy distribution. 

Diagram of repeated convolution of a next scanline error impulse response function 
20 approximating a Cauchv distribution 

The next scanline error impulse response function describes the contribution of 
pixel-only errors of pixels of a scanline to the neighbourhood errors of pixels of the next 
scanline just prior to processing the next scanline. 
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Similarly, the n th next scanline error impulse response function can be defined as 
the contribution of pixel-only errors of pixels of a scanline to the neighbourhood errors of 
pixels of the next scanline just prior to processing the 11 th next scanline. 

Fig. 13 shows the n* next scanline error impulse response function where the 
5 next scanline error impulse response function approximates a Cauchy distribution. It 
shows the contribution of the pixel-only error of a reference pixel to the neighbourhood 
error of a pixel of interest, n scanlines below the reference pixel just prior to processing 
the scanline containing the pixel of interest. Fig. 13 shows a reference pixel 1308 and a 
pixel of interest 1316 separated from the reference pixel 1308 by a distance "k" depicted 
10 by an arrow 1310, in a horizontal direction, and by a distance "n" depicted by an arrow 
13 12 in a vertical direction. 

It is noted that the next scanline 1318 is the scanline immediately following the 
scanline in which the reference pixel 1308 is located. The function 1302 is the 
approximate next scanline impulse response function which determines distribution of 
15 error from the reference pixel 1308 to the next scanline 1318. The function 1304 is the 
approximate second next scanline impulse response function representing distribution of 
error from the reference pixel 1308 to scanline 1320. The function 1306 is the 
approximate n th next scanline impulse response function representing distribution of error 
from the reference pixel to scanline 1322. A bold arrow 1314 represents a radial 
20 separation between the reference pixel 1 308 and the pixel of interest 1316, noting that the 
radial separation is provided by the relationship r=(n 2 '+k 2 ) H , and the approximate impulse 
response is represented by (l/^n/r 2 . 

In this diagram, and hereafter, the Cauchy distribution is referred to with a 
spread parameter, b, equal to 1. Spread parameters close to 1 produce well spread sparse 
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halftone patterns; whereas spread parameters further removed from 1 generate poor 
quality halftone patterns. 

In the diagram, the identical fan of arrows on each scanline is intended to portray 
that the same next scanline impulse response function applies for the processing of each 
pixel of each scanline. 

As has been seen, for the next scanline error impulse response function, h„ CXJ , 
given by (23): 

hnex t [k]* n « (n / 7i) 1 / (n 2 + k 2 ) - (n / tz) 1 / r 2 (1 00) 

So, the cumulative effective of distributing error according to a next scanline 
function approximating (1 / n) 1 / (1 + k 2 ) is that the impulse response of error at the 
reference pixel on the pixel of interest, prior to processing the scanline of the pixel of 
interest, is approximately (n / n) (1 / (n 2 + k 2 ) = (n / n) 1 / r 2 . Here r is the distance 
between the reference pixel and the pixel of interest. 

That is, the pixel-only error at a reference pixel contributes to the neighbourhood 
eiTors of pixels of any succeeding scanline, just prior to processing that scanline, 
approximately in mutual proportions of 1/r 2 . The factor (n / n) simply scales all the 
contributions of the pixels of the particular succeeding scanline. 

It is interesting to note that the impulse response function, c / r 2 , with c a 
constant, is the only radially symmetric function in 2 dimensions for which the impulse 
response on a region is independent of scale. 

To see this, the impulse response on an infinitesimal region, dA, can be 
described as (dA . c/r 2 ). With change of scale by a factor, a, so that horizontal and 
vertical separations x and y instead measure ax and ay, the area of the same infinitesimal 
region now measures a 2 dA, and the impulse response on the infinitesimal region is given 
by: 
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a 2 dA. c / (otr) 2 = dA . c/r 2 (101) 

which is unchanged. 

Graphs of self convolutions of the next scanline e rror impulse response function 

Figs. 14, 15, 16 show graphs of self convolutions of the next scanline error 
impulse response function for: 

1 . Floyd Steinberg error diffusion; 

2. the error diffusion of U.S. Patent No. 5,353,127 (Shiau & Fan); and 

3. where the function is a scaled sampling of a Cauchy distribution. 
Each of these graphs are discussed in relation to the halftone patterns generated 

by corresponding error diffusion algorithms performing monochrome bi-level halftoning 
of 8 bit per pixel image data. 

Fig. 14 shows a graph 1400 of self convolutions of the next scanline error 
impulse response function for Floyd Steinberg error diffusion. The graph in Fig. 14 is 
plotted in regard to an abscissa 1404 which is defined in terms of "horizontal pixel offset" 
from a pixel of interest, and an ordinate 1402 measured in terms of impulse response 
amplitude. A legend 1406 is provided in regard to "next scanline response" (ie. self- 
convolution of degree 1) 1408, "self convolution" (ie. self-convolution of degree 2) 1410, 
"self convolution order 3" 1412, "self convolution order 4" 1414, and "self convolution 
order 5" 1416. In Fig. 14, the self-convolutions are not left-right symmetric and spread 
less rapidly than in Fig. 16 where they spread in proportion to the degree of convolution. 
The sparse halftone patterns of Floyd Steinberg error diffusion exhibit severe worm 
artifacts and strong left-iight asymmetry. 

Fig. 15 shows a graph 500 of self convolutions of the next scanline error impulse 
response function for the patent of Shiau and Fan. The graph in Fig. 15is plotted in 
regard to the abscissa 504 measured in terms of "horizontal pixel offset" from a pixel of 
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interest, and an ordinate 502 measured in terms of impulse response amplitude. A legend 
506 is provided in regard to "next scanline response" 508, "self convolution" 510, "self 
convolution order 3" 512, "self convolution order 4" 514, and "self convolution order 5" 
516, In Fig. 15, the self-convolutions are nearly left-right symmetric, but spread less 
rapidly than in Fig. 16 where they spread in proportion to the degree of convolution. For 
the error diffusion of U.S. Patent No. 5,353,127, worm artifacts are only apparent for very 
sparse halftone patterns. As the halftone patterns become very sparse, the artifacts 
become severe. The worm artifacts show little left-right asymmetry. 

Fig. 16 shows a graph 600 of self convolutions of the next scanline error impulse 
response function for error diffusion as described in relation to Fig. 13. The graph in 
Fig. 16 is plotted in regard to the abscissa 604 measured in terms of "horizontal pixel 
offset" from a pixel of interest, and an ordinate 602 measured in terms of impulse 
response amplitude. A legend 606 is provided in regard to "next scanline response" 608, 
"self convolution" 610, "self convolution order 3" 612, "self convolution order 4" 614, 
and "self convolution order 5" 616. In Fig. 16, the self-convolutions are left-right 
symmetric and spread in proportion to the degree of convolution. In particular, it can be 
seen that the central peak value decreases as 1/n. For error diffusion as described for 
experiment 3 where the next scanline error impulse response function, h„ext, is a sampling 
of the Cauchy distribution (1/ tc) (1 / (1 + x 2 )), normalised so that 2, lw t [i] = 1, and all 
scanlines are processed left to right, there are no worm artifacts, The minority pixels in 
all sparse halftone patterns are well spread. Also, the sparse halftone patterns show very 
little left-right asymmetry. The occasional slight left-right asymmetry can be removed by 
varying the scanline processing direction. 

As previously stated, it is believed that it is desirable that the next scanline error 
impulse response function be left-right symmetric and that it should spread in proportion 
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to the degree of self-convolution. These graphs provide a picture of how existing error 
diffusion methods fall short in this regard. 

A method for generat ing error diffusion masks for which the next scanline error impulse 
response function approximates a Cauchv distribution 

In implementing error diffusion, it is desirable to minimise the processing per 

pixel. 

As stated previously, the processing model of Fig. 11, where error diffusion 
processing per pixel is described in terras of the current and next scanline error impulse 
response functions, is useful for understanding desirable characteristics of error diffusion. 
However, for implementing error diffusion where the next scanline error impulse 
response function approximates a Cauchy distribution, the processing model of Fig. 11 is 
not very suitable. 

This is because the Cauchy distribution reduces slowly away from the peak 
value, so that a large number of non-zero values of the next scanline error impulse 
response function are required to provide a good approximation, with consequently a 
large amount of processing per pixel. 

By contrast, the processing models of Figs. 9 and 10 are better suited to 
implementing Cauchy error diffusion, because a desired next scanline error impulse 
response function can be well approximated using a comparatively small number of mask 
positions. 

A method is now described for designing an error diffusion mask, for which the 
corresponding next scanline error impulse response function approximates a Cauchy 
distribution. 

The difference is measured between the actual next scanline error impulse 
response function corresponding to an error diffusion mask and the desired target next 
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scanline error impulse response function using an objective function, ObjO, equal to the 
sum of the squares of the differences of the function values. 

ObKmaskcurrQ, mask^C) = ^ GV-^afc] " "nexU^M) 2 (102) 
The error diffusion mask weight values described here were generated for a 
target next scanline error impulse response function of the form of equation 23. (Alterna- 
tively, a target next scanline error impulse response function of the form of equation 19 is 
also appropriate). The actual objective function used to determine the error diffusion 
mask weight values is given by: 

Obj (maskc UTT [] s masWO) = 2-n <= k <= m (Iw«»iM ~ (tanhfa) / tt) / (1 + k 2 )) 2 ( 1 03) 
where; 

N is the number of next scanline mask positions to the left of the current pixel 
M = 30. 

It is desired to find mask^] and mask next [] such that 



£, mask^pJ + Sj mask liext [i] = l 
ObjOmask,^ [], mask MX[ []) is a minimum J 



It is noted from (12), that 



(S + h cm ) is a casual function. 



) 
} 



(104) 



(105) 



where 

For hnext to be left-right symmetric as required, mask next should be an anti-causal 

function. That is, mask positions can be omitted on the next scanline ahead of the current 

pixel. It is only necessary to design an error diffusion mask coefficients for pixel offsets: 

with i > 0 on the curent scanline (mask CUfT []), and 1 

with i^Oonthenextscanline(mask lieAt []) J ^° 6 ^ 
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Avoiding use of mask coefficients for pixel offsets on the next scanline with 
i > 0, avoids error distribution processing associated with those pixel offsets, making for 
more efficient error diffusion implementations. 

An approximate solution to this design problem can be solved by optimisation 
methods, including the method of steepest descent and the conjugate gradient method as 
documented, for example, in the textbook: "Applied Numerical Methods for Engineers" 
by R.J. Schilling & S.L. Harris. This textbook also includes C software for performing 
the optimisations. 

The constraint that the sum of the mask coefficients sum to 1 may be effectively 
removed by treating one of the desired error diffusion coefficients as a dependent 
variable. As a first step, to obtaining an error diffusion mask, it is necessary to decide 
how many (or which) mask positions to use on the current and next scanline. That is, an 
error diffusion mask can be designed to meet a pre-determined choice of mask positions 
and consequently a pre-detennined limit of processing per pixel. 

The coefficients for the error diffusion masks of Figs. 4 and 19 were determined 
as described above and using the conjugate gradient method with the aid of software 
supplied with the above textbook. 

Size of error diffusion masks - achieving goals of halftone image quality and imple- 
mentation efficiency 

It has been found that increasing the number of current and next scanline error 
diffusion mask positions, improves the quality of sparse halftone patterns. Although 
judgement of the quality of halftone patterns is somewhat subjective, it has been found 
that, in order to completely remove worm artifacts and poor spreading of sparse halftone 
patterns for the bi-level halftoning of 8 bit per pixel images, it is necessary to use a mask 
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with the width of the error diffusion mask of Fig. 4. The mask of Fig. 4 has 12 positions 
on the current scanline and 13 positions on the next scanline and is called a "12-13 mask". 

The mask of Fig. 4 is used in the description of the error diffusion method of the 
first embodiment. Smaller masks, which are optimised so that the next scanline error 
impulse response function approximates a Cauchy distribution, may also be used - this 
may be desirable to provide a better trade-off between halftone image quality and 
implementation efficiency. 

A family of error diffusion masks for which the next scanline error impulse response 
function approximates a Cauchy distribution 

The most sparse halftone patterns occur for source image regions with grey value 
just above the minimum grey value and just below the maximum grey value. As the grey 
value moves closer to the middle grey value, the halftone patterns become less sparse. 
For most intermediate grey levels, there is no advantage in using error diffusion masks as 
wide as the 12-13 mask. 

It is possible to perform step 4 of the error diffusion processing per pixel, the 
distribution of error, according to a mask dependent on the pixel grey value, This 
approach is used in the error diffusion method of the second embodiment. 

It has been previously observed that the quality of error diffusion halftoning can 
be improved by selectively using a larger error diffusion mask for low and high grey 
value regions. For example, in the paper "Error Diffusion Algorithm with Reduced 
Artifacts", from the Proceedings of IS&T's 45th Annual Conference, May 10-15, 1992, 
New Jersey, the author, R. Eschbach, discusses use of a combination of 2 error diffusion 
masks. 
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In his paper, Eschbach, also observes that it is necessary that the 2 masks should 
have similar structure so that the parts of an image using both masks do not show 
crossover artifacts. 

The requirement that the next scanline error impulse response function should 
5 approximate a particular Cauchy distribution, provides a design goal for generating a 
family of compatible error diffusion masks. 

The second arrangement uses such a family of masks, with a range of widths, to 
achieve similar halftone quality to that achieved by the first arrangement but with 
significantly reduced average processing per pixel. 
10 Description of the First Arran gement 
Bi-level monochrome halftoning 

Fig. 17 is a block diagram 1700, describing the processing performed per pixel 
according to the first arrangement. The processing determines a halftone output value for 
a pixel from the pixel input value and stored error values which are sums of error 
15 distributed to the pixel from previously processed pixels. An input image value g;j for a 
pixel of interest is provided on a line 1702 to an addition process 1704, from which is 
output, on a line 1706, a modified input image value for the pixel (ij), ie. g_modij, this 
being input to an addition process 1708. Thereafter, a table index value table_mdexjj is 
determined, in respect to a table 1726, as the sum of the modified input image value 
20 g_modij, and the least significant bits of the table index value of a previously processed 
pixel which is provided on a line 1718. In the case of left to right scanline processing, the 
table index value is provided by (108) (see below), and in the case of right to left scanline 
processing, the table index value is provided by (109) (see below). It is noted that the 
table index value table_indeXij is computed as the sum of the modified input image value 
25 for the pixel g_modjj on the line 1706 and the least significant bits of the table index 
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value of the previously processed pixel, this having been provided to a latch 1716 and 
then output on the line 1718. The upper bits of the table index value are used via line 
1712 to index the table. The lower bits of the table index value are saved in a latch 1716 
for use in determining the table index value of the next pixel. 

Considering a row 1733 of the table 1726, where the row 1733 has been defined 
by determination of the upper bits of the table index value table_index ij3 then 
commencing at the left hand side of the table 1726, it is seen that respective values CI, 
C2, Cll, C12 are provided on lines 1729, 1731, .... 1732, and 1734 respectively to 
addition processes 1774, 1766, 1750, and also to a latch 1746, respectively. Again 
considering the row 1733 of the table 1726, it is seen that values associated with the right 
hand side of the table 1726, and in particular with columns nO, nl, n2, nl 1 and nl2 are 
provided on lines 1740, 1738, 1736, 1744 and 1742 respectively to a latch 1784, adder 
processes 1790, 1794, 1711, and 1719 respectively. It is noted that an error diffusion 
mask depiction is provided at the top of Fig. 17. The mask refers to a current pixel 720, 
having on its right hand side twelve successive current scanline pixel positions cl to cl2 
(these being collectively designated as reference numeral 1724), and having, on a next 
line, a set of 13 next scanline pixel positions nO to nl2 (collectively designated as 1722). 

A tandem arrangement of successive latches and adder processes 1746, 1750, 
1754, 1762, 1766, 1770, 1774 and 1776 provide, on a line 1778, a total error 
e_euTTjine i%} , of errors distributed directly to pixel (ij) from pixels of the current scanline, 
which is fed to an adder process 1780. In a similar fashion, a tandem arrangement of 
respective latches and adder processes 1784, 1790, 1792, 1794, 1798 1707, 1711, 
1717, 1719 provide, on lines 1723 or 1721, a total of error distributed directly from pixels 
of the previous scanline, ie. e_prev_Iinej+i 2 . j+i (in the case of right to left scanline 
processing), or e_prey_line 1 .i 2l j+i (in the case of left to right scanline processing). The 



191201 



581700US.doc 



values ejprev_linq. 12ij+1 , or e_prev_line i+12i j+1 are directed, by the lines 1723 or 1721 
respectively, to an error line store memory 1725. The memory 1725 outputs, on a line 
1727, the total error e_prevjineij distributed directly to the pixel (ij) from pixels of the 
previous scanline on a line 1727, this being input to the adder process 1780. The adder 
5 process 1780 outputs, on a line 1782, the neighbourhood error value e_nbrij at the pixel, 
this being the sum of errors distributed directly to the pixel, on a line 1782, which is input 
to the adder process 1704. 

In this implementation, the halftone output value on the line 1730 and the errors 
distributed from a pixel are retrieved from a pre-calculated error distribution and result 
10 table 1726. The error values stored in the table have extra precision to store fractional 
error values. A suitable amount of storage per error value is 16 bits, with the most 
significant 8 bits representing the sign of the error together with the integer part of the 
error as a number in the range (-128,127], and with 8 least significant bits representing the 
fractional part of the error. 
15 The same processing cycle is repetitively applied, pixel by pixel across a 

scanline, to process a scanline of pixels. The scanline processing is repeated for each 
scanline to halftone the image. The scanline processing direction for each scanline may 
either be left to right or right to left, and the scanline processing direction may change 
from one scanline to another. Some desirable halftone pattern randomisation is achieved 
20 by varying the scanline processing direction. 

The halftone output value on the line 1730 for a pixel at position (i, j) (being 
column i and scanline j) is determined according to the pixel decision rule defined by the 
following equation. 

if gjj + e_nbrjj > th then r id - = 255; else r;j = 0 (27) 

25 where: 
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gij is the input image value of the pixel; 
e_nbrij is the neighbourhood error value at the pixel, being the sum of errors 
distributed to the pixel; 

th is a threshold value; and 

ry is the halftone output value of the pixel (0 or 255). 

It is convenient, though not necessary for the threshold value to be 127, as this 
then constrains the total error distributed from a pixel to the range (-128,127] and each 
error value stored in the table can be scaled up by 256 and represented as a signed 'short' 
integer in the range (-128*256, 127*256]. With error scaled up by 256 it is also 
necessary to scale up input values by 256. The table index value is then constrained to 
the range (-128,127] x 256 and can be stored in 17 bits. Choosing the upper 9 bits of the 
table index value to index the table produces good halftone output. It is also possible to 
increase the number of bits of the table index used to index the table, in order to allow 
avoiding the use of the lower bits of the table index value altogether; however, this also 
means a consequent increase in the size of the table. 

Determination of a pixel's halftone output value according to the pixel decision 
rule, as well as the determination of the errors to distribute, is achieved by a table lookup. 
The table lookup is made using the most significant bits of a table index value computed 
as follows . 

Firstly, the modified input image value for the pixel (i j), g_modj j, is computed as 

follows: 

g_modij = gy + e_nbrj j ( 1 07) 

Secondly, the table index value, table_index; j, is computed as the sum of the 
modified input image value for the pixel and the least significant bits of the table index 
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value of the previously processed pixel, this being dependent upon the scan direction as 
follows: 

left to right scanline processing 

table_indexij - g_mod io + least significant bits(table_indexi.u) (108) 
right to left scanline p rocessing 

table_indexij = g_modj j + least significant bits(table_index i+1 d ) ( 1 09) 

The error distribution and result table stores the halftone output value and error 
distribution values corresponding to each possible value of the most significant bits of a 
table index value. 

The halftone output value is determined by a thresholding equation similar to 
(27) as follows: 

if table Jndex_upper_bits > th then out = 255; else out = 0 (l 1 0) 

The error distribution amounts are initially determined as 



e _c[iJ = (table_index_upperjrits - out) * mask em . [i] for 1 <= i <= 12 
e _n[i] = (table_index_upper_bits - out) * mask nexl [i] for 0 <= i <= 12 



(111) 



The rounding errors in representing these error distribution amounts as 16 bit 
values are then distributed amongst the error distribution amounts so that 

Si<=j<_ 12 e_c[i] + Z(K=j<=i 2 e_n[i] = table_index_upper_bits - out (112) 
The error resulting by thresholding a pixel is distributed by retrieving error 
distribution amounts from the error distribution and result table. These values are stored 
and summed with previous error distribution amounts and sums of error distribution 
amounts to deliver the error sum values given by formulae, as a function of scan direction 
as follows: 

left to right scanline processing 



e_curr_Iine j+1 j =e_c[lj. . + e_c[2] M j + .. + e_c[12],. n 
e_prevjine i _ l2ij+1 = eja[12).. +e_n[ll] M . +.. + e_n[0] i . I2j 



(113) 
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right to left scanline processing 

e_currjine. u = e_c[l].. + e_c[2] i+lj +.. + e_c[12] 1+lltj 1 

e _prev_line j+12>j+I = e_n[12]. . +e_n[ll] Wj + + e_n[0] 1+]2ij J (114) 

where e_curr_lineij is the total error distributed directly to pixel (ij) from pixels 
of the current scanline, and e_prey_lmeij is the total error distributed directly to pixel (ij) 
5 from pixels of the previous scanline. 

The error value, ejprev_linei. I2 j+i (left to right processing) or e_j>rev_line,. ( - 1 2j +1 
(right to left processing), is written to line store memory from which it will be retrieved as 
part of processing the next scanline. 

The neighbourhood error value for the pixel (ij), e_nbrjj, is derived as the sum 
10 of the error sum value, e_curr_linejj, generated as part of processing the previous pixel of 
the current scanline, and the error sum value, e_prev_line;j ? retrieved from the line store 
memory. 

Description of a Second Arrangement 
Bi-level monoc hrome halftoning 
!5 The second arrangement is also described for the case of halftoning a 

monochrome image. 

This second arrangement uses a family of error diffusion masks, all specially 
designed so that for each mask, the next scanline error impulse response function, 
corresponding to the mask, approximates the same target function; that target function 
20 being a Cauchy distribution which is sampled and normalised. By these means, the 
generated sparse halftone patterns are free of worm artifacts but at the same time the 
processing involved in the first arrangement is reduced. 

Fig. 19 shows an example of a suitable family 1900 of error diffusion masks for 
the second arrangement, The family 1900 of error diffusion masks comprises a first 
25 column 1902 defining a pixel offset in terms of the variables cl to c!2, and n0 to nl2. 
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Error diffusion masks, depicted by columns 1904 to 1926, provide numerical values 
associated with the aforementioned pixel offsets, the masks being designated as mask 1 to 
mask 12 respectively. 

The pixel offsets in the left column of the table correspond to the mask positions 
shown in Fig. 3. Mask 1 of this family has weights for 1 pixel position on the current 
scanline and 2 pixel positions on the next scanline; mask 2 has weights for 2 pixel 
positions on the current scanline and 3 pixel positions on the next scanline; and so on. 
That is, as the mask number increases, the width of the mask increases. Mask 12 is the 
same mask as that described in Figs. 3 and 4. For each mask of Fig. 19, the current line 
mask positions are those closest to the current pixel, and the next line mask positions are 
those closest to the current pixel which are below or behind the current pixel. 

Each mask of the family was prepared as described above, so that its next 
scanline error impulse response function approximates a target next scanline error 
impulse response function given by: 

hncxuan^M = (tanh(7c) / re) / (1 + k 2 )) (l 1 5) 

The target next scanline error impulse response function is derived from the 
Cauchy distribution (1 / it) 1 / (1 + x 2 ) by sampling and normalisation. 

The halftone output value for each pixel is obtained by error diffusion, with the 
error diffusion mask used to determine error distribution amounts for the pixel being 
selected according to the pixel's grey level. 

Fig. 20 shows a table 2000 with a suitable mapping between grey levels and the 
masks of Fig. 19. The table 2000 comprises two columns 2002 and 2004, the first being a 
mask index running from mask no. 1 to mask no. 12, and the second column 2004 
showing correspondence between grey level values and the various mask indices. Thus, 
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for example, mask index no. 6, designated by a reference numeral 2006, is suitable for use 
with grey levels 16 to 18, 237 to 239, 121, and 134. 

For the masks of Fig. 19 and the grey level to mask mapping provided by the 
table of fig. 20, generated halftone patterns do not show crossover artifacts. That is, the 
transition between use of different masks is not apparent in the halftone output. All 
generated sparse halftone patterns are well spread and worm-free; and the appearance of 
the halftone output is very similar to that of the first embodiment. At the same time, the 
average number of addition operations involved in the second embodiment is significantly 
less than that in the first embodiment. 

Fig. 1 8 shows a block diagram for the processing performed per pixel for the 
error diffusion algorithm of the second arrangement. An input image value gy for the 
pixel being considered is provided on a line 1802, this being directed as depicted by lines 
1806 and 1804 respectively to an adder process 1824 and a mask index table 1808. The 
mask index table 1808 provides, in the first instance, a mask index value which is used, as 
depicted by the arrow 1810 to select a particular error distribution table eg. 1814, 1816, 
and 1818, from a family 1820 of error distribution tables. In addition, the mask index 
table 1808 provides two other numbers, these being the number of mask positions on the 
current scanline, and the number of mask positions on the next scanline, for the particular 
mask corresponding to the mask index value. 

The adder process 1824 provides, on a line 1826, a modified input image value 
g_modjj for the pixel (ij), directing this to an adder process 1828. The adder process 
1828 provides, on a line 1830, and thereafter, a line 1832, a table lookup index 
table_index; i j which is determined, as described in relation to Fig. 17, and in accordance 
with (108) and (109). The table index value table_indexij is computed by the adder 
process 1828 as a sum of the modified input image value gmodg on the line 1826 and 
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the least significant bits of the table index value for the previously processed pixel this 
being provided by a line 1834, a latch 1836, and a line 1838 which feeds into the adder 
process 1828. The most significant bits of the table index tablejndexy on the line 1832 
are used to retrieve a halftone output value for the pixel r,j on a line 1822, as well as error 
5 distribution values on the lines 1840 and 1842. The number of current scanline mask 
positions (not shown) and the current scanline error distribution amounts on the line 1840 
are used to update the error sum values in a current line error buffer 1846. 

The number of next scanline mask positions (not shown) and the next scanline 
error distribution amounts on the line 1842 are used to update the error sum values in a 

10 line store error buffer 1848. The current line error buffer 1846 provides an error sum 
value e_currjineij being the total error distributed directly to pixel (ij) from pixels of the 
current scanline on a line 1862 which is directed to an adder process 1860. The line store 
error buffer 1848 provides an error sum value e_prov_linck.i 2 j (for left to right scanline 
processing) on a line segment 1852, or alternatively e_prev_linei +12iJ (for right to left 

15 scanline processing) on a dashed line 1854, for storing in an error line store memory 
1856. The memory 1856 provides, on a line 1858, a value e_prev_line ij5 ie., the total 
error distributed directly to the pixel (i j) from pixels of the previous scanline, which is 
directed to an adder process 1860. The adder process 1860 outputs, a value for the 
neighbourhood error for the pixel (ij), ie., e_nbr u on a line 1844, which is directed to the 

20 adder process 1 824. The nature of the processing performed per pixel is similar to that of 
the first embodiment. The differences in processing between the 2 embodiments are now 
described. 

In the second embodiment, the pixel grey level on the line 1802 is used to 
retrieve a mask index value and 2 other numbers from a mask index table 1808. The 
25 other 2 numbers are the number of mask positions on the current scanline and the number 
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of mask positions on the next scanline, for the mask corresponding to the mask index 
value. The mask index value is used to select an error distribution table from the family 
of error distribution tables 1820. The most significant bits of the table index generated 
for the pixel and the 2 numbers representing the number of mask positions, are used to 
5 retrieve the halftone output value for the pixel and the error distribution values from the 
selected error distribution table. 

The number of current scanline mask positions and the current scanline error 
distribution amounts are used to update the error sum values in a current line error buffer 
O 1846 and to provide the error sum value, e__cujT_line i+ i j (left to right scanline processing) 

nJ io or e_curT_linei.ij (right to left scanline processing) for the next pixel to be processed for 

the scanline. The number of addition operations required for this step is the number of 
current scanline mask positions for the selected mask. This is typically much less than 
the U addition operations required in the first arrangement, 
jjf The number of next scanline mask positions and the next scanline error 

H 15 distribution amounts are used to update the error sum values in a line store error buffer 

1848 and to provide the error sum value, e_prev_linej-, 2si (left to right scanline processing) 
or e_prev_line i+ i2j (right to left scanline processing) for storing in the error line store 
memory 1856. The number of addition operations required for this step is the number of 
next scanline mask positions for the selected mask. This is typically much less than the 
20 12 addition operations required in the first arrangement. 

In summary, the advantage of the second arrangement over the first, is that the 
average number of addition operations required is dramatically reduced. When the 
halftoning method is implemented in software, where each addition operation must be 
performed sequentially, this can significantly improve the execution speed. 
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The described arrangements can be implemented for error diffusion where the 
output halftoned image has more than 2 levels and where the source image is a colour 
image. 

The method of Cauchy error diffusion can be practiced using a general-purpose 
computer system 2100, such as that shown in Fig. 21 wherein the processes of Figs. 9-10, 
17-1 8 may be implemented as software, such as an application program executing within 
the computer system 2100. In particular, the steps of method of Cauchy error diffusion 
are effected by instructions in the software that are carried out by the computer. The 
software may be divided into two separate parts; one part for carrying out the Cauchy 
error diffusion methods; and another part to manage the user interface between the latter 
and the user. The software may be stored in a computer readable medium, including the 
storage devices described below, for example. The software is loaded into the computer 
from the computer readable medium, and then executed by the computer. A computer 
readable medium having such software or computer program recorded on it is a computer 
program product. The use of the computer program product in the computer effects an 
advantageous apparatus for Cauchy error diffusion. 

The computer system 2100 comprises a computer module 2101, input devices 
such as a keyboard 2102 and mouse 2103, output devices including a printer 21 15 and a 
display device 21 14. A Modulator-Demodulator (Modem) transceiver device 2116 is 
used by the computer module 2101 for communicating to and from a coinmunications 
network 2120, for example connectable via a telephone line 2121 or other functional 
medium. The modem 21 16 can be used to obtain access to the Internet, and other 
network systems, such as a Local Area Network (LAN) or a Wide Area Network (WAN). 

The computer module 2101 typically includes at least one processor unit 2105, a 
memory unit 2106, for example formed from semiconductor random access memory 
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(RAM) and read only memory (ROM), input/output (I/O) interfaces including a video 
interface 2107, and an I/O interface 21 13 for the keyboard 2 102 and mouse 2103 and 
optionally a joystick (not illustrated), and an interface 2108 for the modem 2116. A 
storage device 2109 is provided and typically includes a hard disk drive 21 10 and a 
floppy disk drive 2111. A magnetic tape drive (not illustrated) may also be used. A CD- 
ROM drive 2112 is typically provided as a non-volatile source of data. The 
components 2105 to 2113 of the computer module 2101, typically communicate via an 
interconnected bus 2104 and in a manner which results in a conventional mode of 
operation of the computer system 2100 known to those in the relevant art. Examples of 
computers on which the described arrangements can be practised include IBM-PC's and 
compatibles, Sun Sparcstations or alike computer systems evolved therefrom. 

Typically, the application program is resident on the hard disk drive 21 10 and 
read and controlled in its execution by the processor 2105. Intermediate storage of the 
program and any data fetched from the network 2 120 may be accomplished using the 
semiconductor memory 2106, possibly in concert with the hard disk drive 2110. In some 
instances, the application program may be supplied to the user encoded on a CD-ROM or 
floppy disk and read via the corresponding drive 2112 or 2111, or alternatively may be 
read by the user from the network 2120 via the modem device 2116. Still further, the 
software can also be loaded into the computer system 2100 from other computer readable 
medium including magnetic tape, a ROM or integrated circuit, a magneto-optical disk, a 
radio or infra-red transmission channel between the computer module 2101 and another 
device, a computer readable card such as a PCMCIA card, and the Internet and Intranets 
including email transmissions and information recorded on websites and the like. The 
foregoing is merely exemplary of relevant computer readable mediums. Other computer 
readable media may alternately be used. 
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The method of Cauchy error diffusion may also be implemented in dedicated 
hardware such as one or more integrated circuits performing the functions or sub 
functions of Cauchy error diffusion. Such dedicated hardware may include graphic 
processors, digital signal processors, or one or more microprocessors and associated 
5 memories. 

Industrial Applicability 

It is apparent from the above that the arrangements described are applicable to 
the image processing industries. 

The foregoing describes only some embodiments of the present invention, and 
10 modifications and/or changes can be made thereto without departing from the scope and 
spirit of the invention, the embodiments being illustrative and not restrictive. 
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Appendix A 

Abstract 

Although error diffusion is a very popular and effective halftoning method it 
generates halftone patterns with a number of shortcomings. 

One of the major shortcomings is the generation of poor sparse halftone patterns 
where the exceptional pixel values, instead of being well separated, are concentrated in 
wavy lines. These patterns are known as "worm artifacts". 

Various modifications to error diffusion have been devised to minimise or 
eliminate worm artifacts. 

This paper describes a novel view of error diffusion where the error distribution 
from a pixel, resulting from the complete processing of a scanline, approximates a 
standard Cauchy distribution, having the form (1/tt)(1/(14-x 2 )). Error diffusion which 
approximates this error distribution, results in halftone patterns which are well spread at 
all scales; that is, sparse patterns are well spread no matter how sparse they are. 

Analysis and experiments are described which support the desirability of this 
error distribution and provide insight into the error diffusion process. The analysis looks 
at the family of distributions generated by self-convolution from the initial distribution of 
error to the next scanline and the rate at which these distributions spread horizontally. It 
is shown that the family of distributions based on the standard Cauchy distribution is 
radially balanced. 

Error diffusion algorithms which achieve close approximations to this error 
distribution require wide error diffusion masks. However, for a given restriction on the 
size of an error diffusion mask, the principle that the error should be distributed to the 
next scanline as a standard Cauchy distribution provides a useful aid in the design of the 
mask. 
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Error diffusion masks designed according to this principle provide a way of 
implementing error diffusion, which is efficient at reducing worm artifacts, or, if 
sufficient mask positions are used, eliminating them. 

A digital image is typically represented as an intensity value or set of intensity 
values at each of a set of spatial locations ("pixels"). 

Digital halftoning is the transformation from one digital image to another, where 
the output image pixel values are members of a very reduced set of possible values, and 
where the output image appears similar to the original input image. This requires a 
collaborative effect of pixel output values, so that the average intensity in regions of the 
output image approximate the average intensity of those regions in the input image. 
Because the output image thereby renders intermediate intensities or "tones", the process 
is called halftoning. 

Digital halftoning is used to render digital images on devices with limited set of 
output values per pixel. Such devices include inkjet printers (where the printer is 
typically places fixed size dots on a rectangular grid of positions) and some computer 
displays with limited memory available for storing colour values. 

There are 2 main classes of practical digital halftoning algorithms: dither matrix 
halftoning and error diffusion halftoning. Dither matrix halftoning is a point process, 
where each pixel is transformed independently of the surrounding input image values. By 
contrast, error diffusion is a neighbourhood process; where a pixel is transformed taking 
into account the error at previously transformed, neighbouring pixels. Broadly, error 
diffusion involves more computation than dither matrix halftoning but is capable of 
producing more visually pleasing results. 
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Error diffusion 

Error diffusion was invented by Floyd and Steinberg [reference 1], 
In most error diffusion algorithms one pixel is processed at a time. The 
processing performed per-pixel is typically as follows: 
5 1. Add the pixel input value and error transferred to that pixel from previously 

processed pixels to obtain a modulated pixel input value. 

2. Use the modulated pixel input value to decide a halftone output value for the 
pixel, typically by a thresholding operation. 

3. Calculate the combined error at the pixel by taking the difference between the 
10 modulated pixel input value and the pixel halftone output value. 

4. Distribute all the combined error at the pixel to a small set of as yet unprocessed, 
neighbouring pixels according to a set of weights known as an error diffusion 
mask. 

Many extensions and enhancements to error diffusion have been proposed over 
15 the years. The paper "Evolution of error diffusion" by Knox [reference 2] provides an 
overview of enhancements to error diffusion together with developments in the 
understanding of error diffusion. 

Typically, error diffusion, as with many image processing algorithms, is 
performed for one row (or "scanline") of the image at a time, from the top to the bottom 
20 of the image, and pixel by pixel along each scanline. 

This paper will discuss a commonly used subset of such error diffusion 
algorithms where error at a pixel is distributed to pixels of the current scanline and pixels 
of the next scanline only. Also, only error diffusion of monochrome images, with a single 
intensity value per pixel, will be discussed. Finally the halftoning is assumed to be bi- 
25 level; that is, each pixel in the output image may have one of 2 possible values. 
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Extensions of error diffusion to process colour images, with multiple levels of output 
intensity are well known. 

Fig. 22 depicts a common method of performing the processing per-pixel for 
such error diffusion. 

5 The diagram shows an input (greyscale) pixel value, gy, for a current pixel at 

column i and scanline j of the image. A modulated greyscale pixel value, g_modulatedy, 
is derived from the greyscale pixel value by addition of a sum of errors from 
neighbouring pixels, e_nbry. 

By a thresholding operation (indicated as a circled step function), the modulated 
10 greyscale value is transformed to a pixel output value (result) ry. The combined error at 
the pixel, e_combinedy is computed as g^modulatedjj — ry; it is the sum of e_nbry and 
the pixel-only error ey = gy - ry. 

This combined error is distributed in full to neighbouring as yet unprocessed 
pixels on the current and next scanline according to an error diffusion mask, defined by 
15 two sets of weights, represented as discrete 1-dimensional sequences, maskc UIT and 
masknct. 

maskcunEk] is the fraction of combined error, e_combinedy distributed directly to 
the pixel which is horizontally offset by k pixels from the current pixel. As error can only 
be distributed to pixels ahead of the current pixel, maskcun-fk] = 0 for k <= 0. 
20 mask„ ext [k] is the fraction of combined error distributed directly to the pixel on 

the next scanline, offset horizontally by k pixels from the current pixel. 

Portions of combined error, each equal to the product of the combined error with 
one of the weights of masker are added into a current line error buffer, each value in the 
buffer corresponding to a pixel on the current scanline ahead of the current pixel. 
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Similarly, portions of combined error, each equal to the product of the combined 
error with one of the weights of mask„ ext are added into a linestore error buffer, each value 
in the buffer corresponding to a pixel on the next scanline. As part of the per-pixel 
processing, one value of the linestore error buffer will be complete as it will receive no 
further contributions, and it is transferred to the error line store, for use in processing the 
next scanline. 

The error line store holds the error values, e_prev_linej j, each being the sum of 
combined errors distributed directly to a pixel of one scanline fiom pixels of the previous 
scanline. As part of the per-pixel processing, one value is read from the line store and one 
value is written to the line store. Consequently, the line store need only be as long as a 
scanline (hence its name). 

The sum of portions of combined error distributed directly to a pixel from pixels 
of the current scanline, e_curr_lineij, is retrieved from the current line error buffer. The 
sum of portions of combined errors distributed directly to a pixel from neighbouring 
pixels, e_nbrjj and used to modulate the pixel input value is obtained as the addition of 
e_curr_lineij, drawn from current scanline error memory and e_prev_linejj, drawn from 
the error line store. 

Fig. 22 is a block diagram of the processing per-pixel for error diffusion which 
proceeds 1 scanline at a time and which distributes error to the current and next scanlines 
only 

Where the mask is fixed for the entire scanline, the computation of the error 
values, e_curr_linej j and e_prev_linei j for a scanline, can be described by the following 
equations: 

e_curr_linej - (e_linej + e_curr_linej + e_prev_linej) * masker 
e_prev_linej+i = (e_linej + e_curr_linej + e_prev_line,) * mask nex t 
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where: 

e_linej is the sequence of error values, ejj, of scanline j, each of which is the 
pixel-only error at pixel (ij) given by (g ; 0 - - ri d ); 

e_curr_linej is the sequence of error values, e_curr_lineij, of scanline j, each of 
5 which is the sum of error distributed directly to pixel (i j) from pixels of the same 
scanline; 

e_prev_liuej is the sequence of error values, e_prev_line jdJ of scanline j, each of 
which is the sum of error distributed directly to pixel (ij) from pixels of the previous 
scanline; 

10 * refers to the convolution operation between 2 sequences 

(f*g)[i] = 2* f[i-k].g[k] 
An alternative model for error diffusion in terms of the current and next scanline error 
impulse response functions 

Error diffusion processing described by Fig. 22 can also be performed according 
15 to a different model (see Fig. 23). According to the error diffusion processing model of 
Fig. 23, at each pixel, a sum of the pixel-only error and the error distributed to the pixel 
from previous scanlines, is distributed using different error distribution functions. 

To distinguish these error distribution functions from the mask distribution 
functions masker- and mask„ e xi we call them "error impulse response" functions, denoting 
20 them hemr and h B ext- 

Fig. 23 is a block diagram of the processing per-pixel for error diffusion in terms 
of error "impulse response" functions 

hcurr, the "current scanline error impulse response function", is a discrete 1 
dimensional signal representing the pattern of distribution of error at a reference pixel to 
25 pixels of the current scanline resulting from complete processing of the scanline. 
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IWt, the "next scanline error impulse response function", is a discrete 1 
dimensional signal representing the pattern of distribution of error at a reference pixel to 
pixels of the next scanline resulting from complete processing of the scanline. 

To achieve the same output as the original model of error diffusion it is required 

5 that 

e_curr__line, = (e_linej + ejprev_linej) * rw 
e_prev_linej + i = (e_linej + ej>rev_linej) * rw t 

Similarly to raaskc Urr , hc un -[k] = 0 for k<= 0, as pixel processing along a scanline 

is causal. 
10 hcun- can be expressed as 

hc UtT = mask cun . + mask^n-* 2 + maskcun-* 3 .. 

where 

the notation refers to self-convolution (n-1) times, so that 
h* 5 =h; h* 2 = h*h; h* 3 = h*h*h; etc 
15 hnext can be expressed as 

hnext - (5 + hc urT ) * mask next 

where: 

5 is the delta sequence <S[k] =* 1 for k=0; and S[k] = 0 elsewhere) 
The error diffusion processing model of Fig. 23 is more general than that of Fig. 
20 22 - any processing which can be performed according to the model of Fig. 22 can also 
be performed according to the model of Fig. 23; but not vice versa. The model of Fig. 23 
is not suitable for practical implementations of error diffusion - the processing of each 
pixel involves distribution of error to too many pixels. However, the model is useful both 
conceptually and for performing error diffusion experiments as will be discussed. 
25 Error diffusion as a 1 -dimensional process and the "influence function" 
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Although error diffusion is clearly 2-dimensional, error diffusion which proceeds 
one scanline at a time from the top to the bottom of the image, can be viewed as a 1- 
dimensional process. 

The entire current scanline is like a wavefront. Behind the wavefront, processing 
5 is complete and the halftone output is fixed; processing at the wavefront is based on a 
memory of all the processing which has taken place; ahead of the wavefront the halftone 
output is untouched, hi this view, the operation of processing a scanline is considered as 
an atomic operation, a black box. The scanline processing operation takes as input a 
scanline of the input image together with a set of errors from processed scanlines. It 
10 generates as output a scanline of the halftone output image together with an updated set of 
errors from processed scanlines for use in processing the next scanline. For error 
diffusion which distributes error to the current and next scanlines only, the set of errors 
from processed scanlines is 1-dimensional. 

This 1-dimensional view of error diffusion is illustrated in Fig. 24. 
15 Fig. 24 shows a 1-dimensional view of error diffusion where scanline error is 

distributed only to the next scanline 
In Fig. 24, 

g_linej is the current (j'th) scanline of pixel input values (gy), 
r_linej is the corresponding scanline of pixel output values (ry), 
20 The scanline of values written to the error line store can be written as 

e_prev_lineg+i = 2 k (e_prevjinei_kj + e^). iwtk] 
The same calculation can be written in terms of 1-dimensional signals as 
e_prev_linej+i = (e__prev_linej + e_linej) * h n „ t 

Assuming the pattern of error distribution to the next scanline does not change 
25 from one scanline to another, we can write: 
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e_prev_lin©j+i = £ „>=i(e_line j+i . n * h ne xt* n ) 

That is, the values in the error line store can also be viewed as sums of pixel-only 
eiror values of previous scanlines, each convolved with the appropriate member of a 
family of distributions derived from the next scanline distribution function by repeated 
5 self-convolution. 

We can also turn this around and quantify how much the pixel-only error at a 
reference pixel contributes to the error line store values of subsequent scanlines. In a 
sense, this also quantifies the contribution of that pixel-only error to the deterrnination of 
the halftone output result of pixels of subsequent scanlines. The contribution of the pixel- 
10 only error at a pixel to the next scanline is described by rwt; to the scanline after that it is 
described by h n <*. * 2 ; and to the nth next scanline it is described by h next * n . 

I use the term "influence function" to describe the contribution of pixel-only 
error at a pixel to the error line store value of subsequent scanlines. 

This picture of the influence of pixel-only error at a pixel on the determination of 
15 pixels of subsequent scanlines does not attempt to account for the following 2 effects: 

1. the processing of a scanline, typically proceeding pixel by pixel from one end of 

the scanline to the other, being left-to-right or right-to-left accumulates errors 

from pixels of the current scanline together with the error distributed to those 

pixels; 

20 2. the (positive or negative) error at a pixel tends to bias the determination of 

pixels, in its vicinity and processed later, to a value which generates an opposite 
signed error, thereby reducing the effect of the error at the original pixel on other 
pixels processed yet later again. 

Effect 1, inevitably introduces some left or right delay shift of the halftone result 
25 for the scanline with respect to the scanline input. However, this effect can be minimised 
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by choice of the error diffusion mask coefficients for the current scanline. Also, the effect 
can be reduced for an image by processing scanlines in different directions. It is very 
common for error diffusion to be performed with alternate scanlines being processed in 
opposite directions; this is called serpentine processing. 

Floyd-Steinberg error diffusion 

The original Floyd Steinberg error diffusion algorithm processes each scanline 
left to right and distributes error according to the set of weights shown in Fig. 25. 

Fig. 25 shows weights used in the Floyd-Steinberg error diffusion algorithm. 
The influence function for the Floyd-Steinberg error diffusion algorithm is 
shown in Fig. 26 as a surface plot. It is clearly not left-right symmetric. Figs. 31-33 
presents another style of pictorial representation of the influence function. The next 
scanline error distribution function or "next scanline error impulse response function" is 
shown for the first 5 subsequent scanlines, all superimposed on the same 2D graph. 

Fig. 26 shows an influence function for Floyd-Steinberg error diffusion 
algorithm. 

Fig. 28 shows halftone output generated by 3 different error diffusion algorithms. 
In each case the output was generated by processing an 8-bit per pixel monochrome 
image of constant grey value 253 (possible grey values being 0 to 255). 

The influence function of Floyd-Steinberg error diffusion shows a bias to the 
right. With increasing scanline separation from a reference pixel, the influence reduces 
least at an angle of about 20° from vertical. The Floyd-Steinberg error diffusion output in 
Fig. 28 shows both start-up worm artifacts where worm patterns are typically horizontal 
lines and more typical worm artifacts which are oriented sloping up to the right. The 
general slope of these non-start-up worm artifacts is approximately perpendicular to the 
direction of least decay in the Floyd Steinberg error diffusion influence function. 
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Shiau and Fan extended distributio n set error diffusion 

In US Patent 5,353,127, "Method for quantization gray level pixel data with 
extended distribution set" [reference 3], and in the paper, "A set of easily impleraentable 
coefficients in error diffusion with reduced worm artifacts" [reference 4], Shiau and Fan 
describe error diffusion with a mask which extends further back than previous masks. 

The mask is shown in Fig. 26 Error diffusion with this mask greatly reduces the 
presence of worm artifacts compared to Floyd-Steinberg error diffusion. 

However, worm artifacts, although much less visible, do in fact remain, as can be 
seen in the middle halftone output of Fig. 28. 

Fig. 27 shows weights used in error diffusion according to Shiau and Fan (US 
Patent 5,353,127) 

The next scanline error impulse response function for the error diffusion of Shiau 
and Fan is much more left-right symmetric than Floyd-Steinberg error diffusion. The 
worm artifacts in the middle halftone output of Fig. 28 are nearly horizontal. 

Fig. 28 shows error diffusion outputs for a region of greylevel 253, for the three 
error diffusion methods indicated. 

Error diffusion where the next scanline error impulse response approximates a Cauchv 
distribution 

Fig. 29 shows a wide error diffusion mask, for which the next scanline error 
impulse response approximates a Cauchy distribution 

Fig. 30 shows mask weights for the various mask positions in the mask of Fig. 

29. 

The error diffusion mask coefficients of the wide mask of Fig. 29 were obtained 
by an optimisation process to minimise a difference between the actual next scanline error 
impulse response and the sequence (I/71) tanh(7t) l/(l+k 2 ). 
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Deterniining whether worm artifacts are absent or present in a sparse halftone 
partem is somewhat subjective. Nevertheless, it can be reasonably said that using the 
wide error diffusion mask of Fig. 29, any sparse halftone pattern generated for an 8-bit 
per pixel image greylevel is well spread and worm free. 

The halftone output at the right of Fig. 28 was generated with this wide mask and 
shows no worm artifacts. 

Figs. 31-33 show next scanline error impulse response functions for Floyd 
Steinberg, Shiau and Fan, and Cauchy methods. 

Desirable characteristics of the next scanline error impulse response function 

As discussed, error diffusion where error is distributed to the current and next 
scanlines only according to an error diffusion mask, can also be viewed in terms of the 
model of Fig. 23. 

According to the error diffusion model of Fig. 23, the input to the processing of a 
scanline, is actually a scanline of values, each being the sum of an input greyscale value 
and an error Hnestore value. It is reasonable to view the role of the processing of each 
scanline as being to deliver the best halftone results it can, given its input. In many 
successful error diffusion algorithms, the line of input greyscale values is not pre-biased 
to the left or right. It is reasonable to consider that the delay shift effect in processing a 
scanline, as described above, is not particularly significant, So it reasonable to expect that 
just as the current scanline input need not be shifted, so the memory of previous lines 
error, and the next scanline error impulse response function should not be significantly 
biased to left or right. 

So, we say, the next scanline error impulse response function, at least to a first 
approximation, should be left-right symmetric. 
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Also, it is also reasonable to expect that the contribution of error at a reference 
pixel, as measured by the next scanline error impulse response function, to the 
determination of pixel output for pixels on the next scanline, should reduce with 
increasing distance from the reference pixel. Summarising, we say, the next scanline error 
impulse response function should be left-right symmetric and generally bell-shaped. 

With error diffusion, the principle that error should be distributed in full, without 
amplification or attenuation, is almost universally followed. Use of the principle ensures 
that the average intensity of regions of the output image tend to match the average 
intensity of regions of the input image. 

The requirement that error is distributed in full is expressed by the following 
equation concerning error diffusion mask coefficients: 

Sk mask culT [k] + maskncxtfk] = 1 

This requirement places the following restriction on the next scanline error 
impulse response function 
2fchn e xtM = l 

That is, all error must be distributed to the next scanline. Consequently all error 
is distributed to any subsequent scanline as measured by the nth next scanline error 
impulse response function. 

£ k h„ext* n [k] = 1 

The next scanline error impulse response function of the error diffusion of Shiau 
and Fan satisfies the desirable characteristics of being approximately left-right symmetric 
and generally bell-shaped. 

One might suppose that the reason the error diffusion mask of Shiau and Fan 
generates worm artifacts (albeit much diminished compared to Floyd-Steinberg error 
diffusion) is due to the limited width of the mask. 
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Halftone output has been observed which was obtained by error diffusion using 
arbitrarily wide masks of the same form as the Shiau and Fan mask, as specified below. 

Specification of an arbitrarily wide mask of the same form as the Shiau and Pan 

mask: 

N is an arbitrarily large integer, N > 0 
maskc UIT [l] = Vz 

maskeurrlX) =0 for k * 1 

maskn^t-N] - (l/2) !+N 
mask next [-k] = (1/2) 2 "* for N < k <- 0 
mask next [k] = 0 for k elsewhere 

Note that as N increases the next scanline error impulse response function 
corresponding to this mask tends to the discrete, left-right symmetric, next scanline 
response function: 

W[kJ = 1/3 (l/2)« 

In fact, no matter how large N is made such masks still generate worm artifacts 
for sparse patterns corresponding to the most extreme 8-bit greylevels. 

An ideal error distribution for avoiding worm artifacts 

The results of further error diffusion experiments described below provide 
evidence that there is an ideal error distribution for the class of error diffusion algorithms 
which distribute error to the current and next scanlines only, according to a fixed mask. 

This ideal distribution is the member of the family of Cauchy distributions 
known as the standard Cauchy distribution, and is given by 

f([x) = (l/7r)(l/(l+x 2 ) 
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To avoid worm artifacts, not only should a mask be sufficiently wide, but it 
should also provide a next scanline error response function approximating the standard 
Cauchy distribution. 

The Cauchv distribution 
5 The family of Cauchy distributions is given by 

f(x) = (l/7t)(n/(n 2 -Kx-t) 2 ) 
where t is a location parameter and n is a scale parameter. For each distribution of the 
family, the area under the curve is 1. From this point on only members of the family of 
Cauchy distributions will be discussed which are centrally located (t = 0) and which have 
10 integer values of n. 

Samples of members of this family of Cauchy distribution (with n = 1,2,3,4,5) 
are shown in Fig. 33. The Cauchy distribution itself is bell-shaped, extends without limit 
along each end of the x axis, and has significant weight in the tails - its absolute moment 
(expected absolute value) is not finite as the corresponding integral diverges. 
15 The Cauchy distribution also possesses the remarkable property that its self- 

convolutions are also Cauchy distributions. Indeed, 

[(1 / n) (1 / (1 + x 2 ))] * n = (1 / n) (n / (n 2 + x 2 )) [derived in Appendix B] 
Consequently, a next scanline error impulse response function approximating the 
standard Cauchy distribution, has its n'th scanline error impulse response function 
20 approximating the central Cauchy distribution with scale parameter n. 

Members of the family of Cauchy distributions are copies of each other, obtained 
by inversely scaling in the 2 axes. In particular, the central Cauchy distribution with scale 
parameter n is obtained from the standard Cauchy distribution by stretching the x axis by 
n and by scaling the ordinate axis by multiplying by 1/n. The total area under the curve is 
25 unchanged by this double scaling. 
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(1/n) * 1 / (1 + (x/n) 2 ) = n / (n 2 + x 2 ) 

So it is seen that, the members of the family of central Cauchy distributions with 
n = 1, 2, 3, spread horizontally in proportion to the scale parameter n. 
Experiments exploring the role of the shape of the next scanline distribution function in 
5 generation of worm artifacts 

A range of experiments has been conducted, where the shape of the next scanline 
error impulse response function was varied and the halftone output was observed. 

The error diffusion processing model of Fig. 23 was used to set guidelines and to 
specify the error distribution for these experiments. 
XO As noted above, desirable characteristics for the next scanline function are that 

1. Z k lw[kJ = l 

2. hiexi is left-right symmetric and monotonically decreases from a central peak 
value at hn^fO] 

A further guideline was set that the current scanline error impulse response 
15 function should have a similar shape to one side of the next scanline error impulse 
response function, being obtained from it by multiplying by a step function. 

The actual method of performing the error diffusion processing was according to 
the model of Fig. 22for experiment 1 (as the corresponding mask is straightforward to 
calculate in this case) and was according to the model of Fig. 23 for the remaining 
20 experiments. 

When the model of Fig. 23 was used, the current and next scanline error impulse 
response functions were truncated by setting coefficients beyond a certain horizontal 
offset to zero so as to keep the required processing time manageable. The number of non- 
zero coefficients was still kept high to avoid misleading results. 
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The next scanline error impulse response function is a sampling of the 2-sided 
exponential distribution: 

hc U rr[k] = a . w* for k >= 1 

b-ncxtM = ((1 - W) / (1 + W» . W 1 *' 

where a is a positive constant and 0 < w < 1 

When a = 1 and w = 54, this form of error distribution matches that already 
mentioned as an extension of the error diffusion of Shiau and Fan. In that case, halftone 
output at extreme grey levels suffers from worm artifacts, where the horizontal separation 
of minority pixels is too small. 

These worm artifacts can be reduced by increasing either of the parameters a or 
w. However, modifying these parameters so that minority pixels at the extreme grey 
levels are better separated, introduces other unpleasant artifacts, including artifacts where 
for other grey levels the vertical separation between minority pixels becomes too small. 

Experiment 2 

The next scanline error impulse response function is a sampling of the gaussian 
distribution; 

IWk] = a . exp(-7t(k/b) 2 ) for k >= 1 

hnextM = C . exp(-71 (k/b) 2 ) 

where a, b, c are positive constants, with c chosen so that E k h ncxt [k] - 1; 
b is the standard deviation of the corresponding continuous distribution 
When a = '/ 2 and b = 4, halftone output at extreme grey levels suffer from worm 
artifacts. These artifacts can be reduced by increasing a or b. However, modifying these 
parameters so that minority pixels at these grey levels are better separate, introduces other 
unpleasant artifacts, including where for other grey levels the vertical separation between 
minority pixels becomes too small. 
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Experiment 3 

The next scanline error impulse response function is a sampling of the Cauchy 
distribution: 

hcurrM = a . (1 / (b 2 + k 2 )) for k >= 1 
5 rwt[k] = c . (1 / (b 2 + k 2 )) 

where a, b, c are positive constants, with c chosen so that 2 k h ne xr[k] = \; 
b controls the degree of spread of the distribution 

When a = 1 and b = 1, the halftone patterns for all grey levels with sparse 
halftone patterns are well spread, showing no worm artifacts. Also the general quality 
10 (absence of unwanted textures) of the halftone output is high across all grey levels. 

For 8-bit per pixel monochrome images, the extreme sparse halftone patterns 
correspond to grey levels 1 and 254. By allowing fractional grey levels less than 1 the 
behaviour of error diffusion can be observed for halftone patterns which are much more 
sparse than the halftone patterns for grey level 1 . 
15 It is noteworthy that the very sparse halftone patterns for sub- unity fractional 

grey values are also well spread when the Cauchy distribution is used for the current and 
next scanline response functions. That is, error diffusion processing using the Cauchy 
distribution appears capable of generating well spread halftone patterns no matter how 
sparse the patterns. 

20 The satisfactory output provided in this case also provides some vindication for 

the choice of a current scanline error impulse response function as a copy of the next 
scanline error impulse response function obtained by multiplication by a step function. 

When b is increased well above 1 the vertical separation of minority pixels in 
sparse halftone patterns is too small; when b is decreased well below 1 the horizontal 

25 separation of minority pixels in sparse halftone patterns is too small. 



191201 



581700US.doc 



-92- 

Experiment 4 

The next scanline error impulse response function is a generalisation of the 
Cauchy distribution: 

hcutrM = a . 1 / (b 2 + kY for k >= 1 
lWk] = c. 1 /(b 2 + ]^f 

where a, b, c, p are positive constants, with c chosen so that£ k IwM = 1 

With p set above 1, the observations are similar in character to those for 
experiments 1 and 2. 

With p set below 1 (but necessarily greater than 1/2), the halftone patterns for 
very low or very high grey levels are well spread horizontally, but the horizontal 
separation for intermediate grey values is too small. 
Conclusions and interpretations from the exp eriments 

Of the various next scanline error impulse response functions considered in the 
above experiments, the only form of the next scanline error impulse response function 
which delivers sparse halftone patterns which are well spread for all grey levels is the one 
which approximates the standard Cauchy distribution. 

A next scanline error impulse response function approximating the standard 
Cauchy distribution successfully shapes the error diffusion influence function to achieve a 
suitable balance between providing sufficient horizontal and vertical separation between 
minority pixels, for all grey levels. 

The poorly separated sparse halftone patterns produced by other next scanline 
error impulse response functions can be explained in general terms by considering the 
shape of the corresponding influence function, as follows. 

When the influence function is not horizontally spread sufficiently (too 
concentrated centrally), the disincentive to place minority pixels reduces too rapidly with 
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horizontal separation, resulting in too little horizontal separation between minority pixels, 
causing worm artifacts. When the influence function is spread too much, the disincentive 
to place minority pixels does not reduce rapidly enough with horizontal separation 
causing minority pixels to be placed with a horizontal separation which is too wide and 
consequently with a vertical separation which is too small. 

As the input image grey level approaches the extreme grey levels and halftone 
patterns become more sparse, so the shape of the influence function at more distant 
scanlines becomes important. When the next scanline error impulse function spreads with 
self-convolution at a rate less than in proportion to scanline separation (experiments 1,2, 
and experiment 4 with p above 1), it is possible for the influence function to be too spread 
for near scanlines and not spread enough for distant scanlines. When the next scanline 
error impulse function spreads with self-convolution at a rate greater than in proportion to 
scanline separation (experiment 4 with p below 1), it is possible for the influence function 
to be not spread enough for near scanlines and to be too spread for distant scanlines 
Desirability that the next scanline error impulse response function should spread in 
proportion to the degree of self-convolution 

As previously noted, the Cauchy distribution spreads with self-convolution in 
proportion to the degree of self-convolution. By contrast, a gaussian distribution spreads 
in proportion to the square root of the degree of self-convolution. 

It is proposed that a necessary requirement for well spread sparse halftone 
patterns for error diffusion, where error is distributed to the current and next scanlines 
only and the pattern of distribution is fixed, is that the next scanline error impulse 
response function should spread in proportion to the degree of self-convolution. 

In fact, the requirement that the next scanline error impulse response funchon 
should spread in proportion to the degree of self-convolution, when considered with the 
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other desirable characteristics listed previously determines that the next scanline error 
impulse response function must be a Cauchy distribution. 

Specifically, a positive symmetric continuous impulse response, which, 
preserves the magnitude of the impulse, spreads in proportion to the degree of self- 
convolution and which acts as a smooth (continuous) low pass filter, dampening all non- 
zero frequencies must be a Cauchy distribution [see Appendix B], 

Another view of the Cauchy distribu tion and the desirability that the next scanline error 
im pulse response function should be radia lly balanced 

Consider the 2-dimensional field of vectors oriented away from a point source 
and of magnitude V(tzt), where r is the distance from the point source. This vector field is 
depicted in Fig. 34. A line, Li. at distance 1 from the point source, O, is shown, defining 
an x axis. A parallel line, L„ 5 is shown at distance n from O. 

Fig. 34 shows the standard Cauchy distribution and the l/(wr) vector field. 

The flux of a 2-dimensional vector field across a line is defined as the line 
integral of the component of the vector field perpendicular to the line. The flux of the 
vector field across line segment s-M: of line Lj and across line segment ns-»nt of line 1^ 
are given directly above (see Appendix B for evaluation of the integral) 

That is, the flux density function along the line L] is the standard Cauchy 
distribution, and the flux density function along the line Ln is the central Cauchy 
distribution of scale n. 
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The flux of the vector field across the line segment s-* t of Li and the line 
segment ns-mt of U are the same, being the flux within the radial segment defined by 
the rays 0-*s and 0-»t and subtending angle e«. Moreover, the flux across any other 
segment of Li, u-»v which subtends the same angle, matches the flux across the segment 
s-M. In this sense, the standard Cauchy distribution is "radially balanced". Indeed, 
another definition of a Cauchy distribution is that probability distribution of horizontal 
distances at which a line, through a fixed point, tilted at a random angle cuts the x-axis. 

So it is seen that, the family of Cauchy distributions (1/n) n/(n 2 +x 2 ) describes the 
flux density of the I/tit) vector field across the parallel lines L„. Consequently, when 
error is distributed to the next scanline according to the radially balanced, standard 
Cauchy distribution, and to subsequent scanlines by repeated self-convolution, it is seen 
that the strength of the distribution within any radial segment is maintained from one 
scanline to the next. Neither does the strength of the distribution concentrate towards the 
centre; neither does it drift to the extremities. 

Design of error diffusio n masks bv optimisation of a cost function. 

Our goal is to design error diffusion masks for which the corresponding next 
scanline error impulse response function approximates the standard Cauchy distribution. 
To achieve this it is necessary to deal with discrete sequences rather than continuous 
distributions. 

The sequence obtained by sampling the standard Cauchy distribution and 
multiplying by the constant tanh(Tt), has its sum of coefficients equal to 1 . 

S k (tanh(Tt) / 7t) 1/(1 + k 2 ) = 1 [see Appendix B] 

The sequence 

h[k] = (1 + (-l) k -V*) / (1 + k 2 ) 

satisfies 
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h[k] * n = (n/ic) (1 + (-ly-V™) / (n 2 + k 2 ) [see Appendix B] 
I k h[k] = 1 [see Appendix B] 

Now, e* * .0432 is close to zero, and tanh(u) * .9963 is very close to 1. So, 
setting a target next scanline error impulse response function equal to 
5 iw^tM - (tanhfr) / «) l/d + k 2 ) 

means that 

2k fincxtjarget M = 1 

That is, the sequence h„««_u.rg« is a sampling of the standard Cauchy distribution 
10 multiplied by tanh(7i) and its self-convolutions are approximately samplings of the self- 
convolutions of the standard Cauchy distribution. 

Having determined h WBet as a suitable next scanline error impulse response 
function, it is then required to determine a corresponding error diffusion mask. 
We refer back to the equation 
15 Iw = (5 + lw) * mask ncKt 

and note that (8 + hc UIT ) is a causal sequence. For the next scanline error impulse 
response function, h nex t, to be left-right symmetric, it is necessary that the sequence 
maskncxt be a scaled anti-causal copy of (5 + h clirr )- This means, that the error diffusion 
mask coefficients, mask^M, should only be non-zero for k < 0. 
20 For an error diffusion mask of a specified width, we seek mask coefficients 

which minimise a difference between the actual next scanline error impulse response 
function corresponding to those mask coefficients and the target next scanline error 
impulse response function. The width of the error diffusion mask is described by 
maximum offset non-zero mask positions for the current scanline, N c and next scanline, 
25 N n , giving 
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maskcuirfk] * 0 for 1 < k <, N c and 
mask ncxt [kj * 0 for -N„ < k < 0. 



We seek (N c + N n + 1) mask coefficients subject to the requirement that 
5 Sk maskcurrM + 5* mask ncxt [k] = 1 

which minimise a suitable cost function such as 

CoStCmaskcuwt], mask next []) = 2 NS M 0W_.ch«l[k] - h nBXlJare ei[k]) 2 

where 

M specifies a suitable large bound for sequence offsets, and as above, 
10 h. n ext_actuai = (8 + maskc UIT + mask Curf * 2 + maskcuir* 2 + mask cwr * 3 ..) * mask^ 

hnexLlareeitk] = (tanh(Tt) / 7C) 1/(1 + k 2 ) 

An approximate solution (mask coefficients which deliver a cost difference near 
the minimum) can be obtained using optimisation methods including the method of 
steepest descent and the conjugate gradient method [see reference 5]. The coefficients for 
15 the error diffusion mask of Fig. 29, were obtained using the conjugate gradient method 
and C software supplied with that textbook. 

A family of error diffusion masks and more efficient error diffusion algorithms 

Using the above method, error diffusion masks have been prepared for a range of 
widths as shown in Fig. 35, Mask 1 has 3 mask positions, one on the current scanlme and 
20 2 on the next scanlme, mask 2 has 2 mask positions on the current scanline and 3 on the 
next scanline; and so on. The mask pixel offsets are labelled cl to cl2 and nO to nl2 and 
correspond to the mask positions shown in Fig. 29. Mask 12 of Fig. 35 is the same mask, 
pixel offsets and mask coefficients, as the mask in Fig. 29. 
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Fig. 35 shows error diffusion masks for a range of widths, optimised so that the 
next scanline error impulse response function approximates the standard Cauchy 
distribution 

Many error diffusion masks have been proposed which distribute error to several 
5 scanlines. For each extra scanline in the error diffusion mask, an extra error line store is 
required to implement the algorithm. While the error diffusion discussed in this paper, 
has the advantage that it combats worm artifacts without a cost in additional line stores, 
an error diffusion mask such as that of Fig. 28 is very wide, much wider than error 
diffusion masks used in practice. The processing performed per pixel relates directly to 

10 the number of error diffusion mask positions. The large width of the mask of Fig. 29 
would prevent its use in many applications. 

For bi-level halftoning, the most sparse halftone patterns occur for input image 
regions with grey levels close to the rninimum or maximum intensity. As the grey level 
moves closer to the centre grey level, halftone patterns become less sparse. 

15 It has previously been observed that the equality of error diffusion halftoning can 

be improved by using a larger error diffusion mask for the grey levels with more sparse 
halftone patterns. For example, in the paper "Error diffusion algorithm with reduced 
artifacts" [reference 6], Eschbach, discusses using a small or a large error diffusion mask, 
depending on the input pixel grey level. 

20 Similarly, for most intermediate grey levels, it is found that there is no 

significant advantage in the quality of halftone output obtained by using the wide error 
diffusion masks of Fig. 35. For most grey levels, masks of Fig. 35 with only a small 
width, can be used without any noticeable artifacts. 

In his paper, Eschbach also noted that the 2 masks should have similar structure 

25 to minimise crossover artifacts for image regions where both masks are used. 
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In a similar vein, it is possible to select the error diffusion mask, based on the 
input image grey level, from a family of masks of varying width, all designed according 
to the principle of minimising the difference between the corresponding next scanline 
error impulse response function and the standard Cauchy distribution. The table of Fig. 
36 defines a mapping between grey levels and the error diffusion masks of Fig. 35. 

By performing error diffusion, where the mask is selected according to the input 
pixel grey level and the table of Fig. 36, very similar quality halftone output is achieved 
to that achieved by use of the mask of Fig,. 29 for all input pixel grey levels. 

In this way, a considerable saving of average processing peT pixel can be 
achieved. 

Fig. 36 shows mapping between grey levels and masks of Fig. 35 

Alternatively, simply using a moderately wide mask from the table of Fig. 35, 
provides good reduction in worm artifacts for a modest increase in processing. 
Conclusio n and final notes 

A principle has been proposed for the design of error diffusion masks which 
distribute error to the current and next scanlines only. The principle is that the error 
distribution to the next scanline resulting from the complete processing of a scanline 
should approximate the standard Cauchy distribution. 

By careful design of the error diffusion mask according to this principle, and 
with no other change to the error diffusion algorithm, the halftone output is free of worm 
artifacts, no matter how sparse the halftone patterns. 

Experimental results have been described which support the principle. Analysis 
has been presented as to why such a pattern of error distribution is desirable. The analysis 
is based on considering the family of distributions generated from the next scanline error 
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distribution by self-convolution. The family of distributions generated from the standard 
Cauchy distribution is radially balanced. 

The names Cauchy error diffusion or radially balanced error diffusion are 
suggested for error diffusion where the mask is designed according to this principle. 
5 Actual masks of varying width designed according to the principle are reported. 

While wide masks are required to eliminate worm artifacts, less wide masks or selective 
use of masks of varying width provide implementation options with less processing. 
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Appendix B 

[(1 / w) (1 / 0 + x 2 ))] * n = (1 / n) (n / (n 2 + x 2 )) 

A positive symmetric continuous impulse response, which, preserves the 
magnitude of the impulse, spreads in proportion to the degree of self-convolution and 
which acts as a smooth (continuous) low pass filter, dampening all non-zero frequencies 
must be a Cauchy distribution. 

Integral s->t of Cauchy distribution is 8st/w 

E k (tanh(7t)/7T) 1/(1 + k 2 ) = l 

For, 

h[k] = (l/7i)(l+(-l) k - 1 e' l )/(l+lc 2 ) 
h[k] *" - (n/«) (1 + (-l^'O / (n 2 + k 2 ) 
2 k h[k] = l 
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